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I.   INTRODUCTION 

In  1978,  the  Navy  Personnel  Research  and  Development 
Center,  located  at  San  Diego,  California,  formulated  "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload"  in  the  11th  Naval  District.   At  the  time  of  its 
formulation,  this  model  was  designed  to  forecast  the  workload 
of  twenty  shore-based  activities  within  the  11th  Naval  Dis- 
trict.  This  thesis  is  the  result  of  a  parametric  sensitivity 
analysis  of  that  model.   In  the  sensitivity  analysis,  the 
author  examined  the  output  levels  of  the  model  when  the  input- 
output  coefficient  matrix  was  changed  by  a  small  amount  in 
each  row. 

Input-output  models  may  be  categorized  as  being  descrip- 
tive or  prescriptive  in  nature  depending  upon  1)  the  goal  or 
purpose  of  the  model,  2)  the  applied  use  of  the  model,  3)  the 
structure  of  the  model,  and  4)  other  related  items.   Within 
this  broad  categorization,  it  is  useful  to  classify  models  in 
terms  of  their  use.   Input-output  models  may  be  classified  as 
models  of  pure  logical  consistency  or  models  of  pure  predic- 
tion.   A  model  of  pure  logical  consistency  deals  only  with 
the  model's  logical  consistency.   It  addresses  the  question 
of  whether  or  not  certain  events  logically  lead  to  certain 


Jones,  C.R.,  A  Taxonomy  for  Naval  Force  Level  and  Structure 
Models,  working  papers  Naval  Postgraduate  School,  Department 
of  Operations  Research  and  Administrative  Sciences,  Monterey, 
California,  April,  1971. 
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other  events.   There  are  no  empirical  tests  required  to  check 
this  class  of  model.   The  pure  prediction  model  generates 
forecasts  or  predictions  for  a  subsequent  time  period. 
Empirical  verification  of  this  class  of  models  may  be  done  by 
comparing  forecast  and  realized  observations,   "A  Regional 
Input-Output  Model  For  Forecasting  Shore-Based  Navy  Workload" 
in  the  11th  Naval  District  is  a  member  of  the  pure  prediction 
class  of  models. 

Section  II  of  this  thesis  provides  insight  into  how  an 
input-output  model  is  devised.   In  particular,  it  shows  how 
the  input-output  coefficients  are  computed  for  each  sector/ 
activity  of  a  general  input-output  model  and  gives  a  specific 
example  of  the  formulation  and  workings  of  a  very  simple  model 
The  mathematics  of  the  general  input-output  model  is  gone 
through  in  detail  and  the  model  is  reduced  in  form  to  a  simple 
expression  in  matrix  and  vector  notation. 

Presented  in  Section  III  is,  "A  Regional  Input-Output 
Model  For  Forecasting  Shore-Based  Navy  Workload,"  the  model 
formulated  by  the  Navy  Personnel  Research  and  Development 
Center.   A  step-by-step  breakdown  of  the  model's  formulation 
and  background  on  how  the  data  were  gathered  are  given.   Many 
years  of  data  collection  and  analysis  went  into  the  formula- 
tion of  the  model.   This  section  presents  the  method  for  aggre- 
gation of  sectors  and  coefficients. 

The  methodology,  chosen  by  this  author,  to  do  the  sensi- 
tivity analysis  on  "A  Regional  Input-Output  Model  For  Fore- 
casting Shore-Based  Navy  Workload"  in  the  11th  Naval  District 
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is  presented  in  Section  IV.   The  decision  to  do  a  parametric 
sensitivity  analysis  on  the  Input-Output  Coefficient  Matrix, 
A,  and  to  only  parametrize  rows  as  opposed  to  rows  and  columns 
or  some  combination  of  rows,  columns,  and  individual  elements 
of  the  input-output  coefficient  matrix  is  explained.   Also, 
the  choice  of  change  values,  for  parametrization,  of  between 
minus  25  percent  and  plus  25  percent  is  explained.   The  pro- 
cedure for  computing  the  output  level  of  the  parametrized 
row/sector  is  outlined  as  is  the  procedure  for  computing  the 
Baseline  Output  Level  used  in  the  calculation  of  percentage 
changes  in  the  output  levels  of  the  twenty  sectors  of  the 
model.   The  crux  of  this  thesis  is  the  sensitivity  analysis 
of  the  output  level  of  a  sector  based  on  a  parametric  change 
to  that  sector's  row  of  the  input-output  coefficient  matrix 
and  the  manner  in  which  this  was  done  is  outlined  in  this 
section. 

The  results  of  the  methodology  outlined  in  Section  IV, 
when  applied  to  "A  Regional  Input-Output  Model  For  Forecasting 
Shore-Based  Navy  Workload,"  are  given  in  Section  V.   Upper  and 
lower  bounds  on  the  percentage  change  in  the  output  level  of 
a  parametrized  sector  are  computed  and  graphed  for  the  twelve 
sectors  this  author  felt  to  be  of  concern.   These  graphs  are 
presented  as  Figures  3  through  14  of  Section  V.   The  lower 
bounds  showed  a  linear  relationship  with  the  delta  values 
used  for  parametrization  to  compute  the  percentage  change  in 
the  output  level  for  a  particular  sector.   The  same  was  true 
for  the  computed  upper  bounds . 
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Conclusions  and  recommendations  are  given  in  Section  VI . 
The  fact  that  the  Public  Works  Center  was  the  most  sensitive 
sector  of  the  model  and  how  this  was  determined  are  presented, 
Inter-sectoral  relationships  are  discussed  and  the  recommen- 
dation to  reevaluate  the  input-output  coefficient  matrix  in 
terms  of  the  zero  entries  is  made. 
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II.   THE  GENERAL  INPUT-OUTPUT  MODEL  AND  RELATED  MATHEMATICS 

The  classical  Leontief  Input-Output  Model  consists  of 
three  basic  elements.   A  square  matrix,  B,  is  used  to  repre- 
sent the  industrial  classification  of  an  economy.   Each  row 
of  the  matrix,  B,  contains  the  amount  of  a  particular  indus- 
try's final  product  that  is  used  as  input  for  producing  the 
final  product  of  industries  represented  in  the  columns  of 
the  matrix.   A  separate  row  vector  is  used  to  denote  the 
amount  of  each  primary  factor,  or  input  such  as  labor,  re- 
quired by  each  industry  in  the  economy.   Finally,  a  separate 
column  vector  is  used  to  provide  the  amount  of  exogenous 
demand  that  exists  for  each  final  product  [11;13]. 

As  a  simplified  example,  consider  the  two-industry  economy 
in  which  industry  one  produces  water  and  industry  two  produces 
electricity.   If  b. .  denotes  the  quantity  of  the  i    industry's 
final  product  used  by  the  j    industry,  then  it  follows  that 
b, -,  represents  the  amount  of  water  used  in  the  production  of 
water  and  b, ^  represents  the  amount  of  water  used  in  the  pro- 
duction of  electricity.   In  a  physical  sense,  b,,  represents 
the  water  used  to  cool  and  supply  evaporators,  other  plant 
devices,  and  other  related  activities  of  the  two  producing 
industries.   In  a  likewise  manner,  b^-,  and  b2  2  represent  the 
amount  of  electricity  used  to  operate  lights  and  other  elec- 
trical components  in  the  water  and  electrical  industries 
respectively.   The  components  b, ,   ^lo'  ^21  '  ^^^  ^22  '^o^prise 
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the  square  matrix,  B.   Inasmuch  as  this  is  a  two-industry 
economy  example,  the  amount  of  labor  required  can  be  repre- 
sented by  Xq,  and  x-2  ^^^   the  exogenous  demand  for  each  final 
product  by  d^  and  d_.   The  total  real  output  of  the  two 
industries  can  be  represented  as  x,  and  x^ .   Figure  1  is 
that  of  the  general  model  for  an  n-industry  economy. 

The  details  of  solving  the  two-industry  economy  problem 
are  presented  as  an  introduction  to  the  methodology  used  in 
working  with  an  input-output  model.   The  solution  procedure 
deals  with  the  problem  of  determining  the  appropriate  levels 
of  output  at  which  the  two  industries  should  operate  such 
that  production  just  satisfies  the  total  demand  for  a  final 
product.   Simply  stated,  equilibrium  exists  in  the  two-industry 
economy  when  the  total  outputs  of  industry  one  and  those  of 
industry  two  are  in  balance  so  that  just  enough  of  each  is 
produced  to  satisfy  both  the  final  demand  and  the  input  re- 
quirements for  each  product.   In  the  general  case,  the  two- 
industry  economy  becomes  an  n-industry  economy  and  the  solution 
is  analogous  to  the  two-industry  discussion. 

To  solve  the  two-industry  problem,  let  x,  and  x^  represent 
the  total  outputs  of  industries  one  and  two  respectively. 
These  quantities  are  represented  as  row  sums  in  Figure  1.   It 
is  assumed  that  the  input-output  coefficients  are  constant 
[10]  and  that  the  transformation  of  inputs  into  outputs  is 
accomplished  in  a  linear  fashion.   From  the  discussion  of  the 
two-industry  economy  it  follows  that  for  each  unit  of  the  j 
commodity  produced,  a  fixed  input  of  the  i    commodity  is 
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Figure  1.  The  Classical  Leontief  Model* 
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required.   In  other  words,  a  unit  of  output  of  water  or 
electricity  requires  a  fixed  input  of  water  or  electricity 
respectively. 

The  input-output  matrix.  A,  may  be  derived  from  the  inter- 
industry relationships.   In  general,  the  input-output  coeffi- 
cients are  defined  by 


a.  .   =  b.  ./x  . , 


where  a. .  represents  the  amount  of  output  of  industry  i  used 
in  the  production  of  one  unit  of  output  of  industry  j,  b. . 
represents  the  input  by  industry  j  to  industry  i,  and  x. 
represents  the  total  real  output  by  industry  j .   Thus ,  the 
input-output  coefficients  for  the  two-industry  problem  are 

^11  =  ^ll^^'l'  ^12  =  ^12/''2'  ^21  =  ^21'^''l'  ^""^  ^22  =  ^22/''2* 
The  total  output  requirements,  x,  and  x^   niay  now  be 

expressed  as: 


^1   =   ^11^1  -^  ^12^2  -^  ^1 


and 


^2   "   ^21^1  "^  ^22^2  "^  ^2* 


Solving  these  two  equations  for  d.  yields 


^1   =   (l-aii)x^  -  a^2^2 
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and 


^2   =   -^21^1  "•  (1-^22)^2- 


Expressing  these  algebraic  expressions  in  matrix  notation  is 
as  follows: 


(1  -  a^^)  -a^2 


-a 


21 


(l-a22) 


X2 

= 

^1 
<^2 

The  above  matrix  notation  may  be  compactly  written  as 

(I-A)x  =  d. 

This  compactly  written  matrix  notation  will  be  used  through- 
out the  thesis  with  A  denoting  the  matrix  of  a. .s,  x  denoting 
the  vector  of  outputs,  d  denoting  the  vector  of  exogenous 
demands,  and  I  denoting  the  identity  matrix  (the  nxn  matrix 
with  Is  on  the  principal  diagonal  and  Os  elsewhere) .   Solving 
the  above  expression  for  x  completes  the  problem  and  it 
follows  that  X  =  (I-A)   d  .   Now,  it  can  be  seen  that  the 
appropriate  level  of  output  for  each  of  the  two  industries 
is  uniquely  determined  for  a  given  exogenous  demand. 

The  assumptions  relating  to  the  Classical  Leontief  Input- 
Output  Model  will  now  be  addressed.   The  first  assumption 
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relating  to  a  single  output  for  each  sector  poses  no  concep- 
tual problems.   However ,  the  assumption  relating  to  the  con- 
stancy of  input-output  coefficients  has  been  the  subject  of 
major  criticism  of  the  model,   A  very  enlightening  discussion 
relating  to  the  constancy  of  the  technological  coefficients  is 
provided  by  Hatanaka  [10] ,   Hatanaka's  first  area  of  criticism 
relates  to  weaknesses  resulting  from  ignoring  certain  factors 
considered  in  production  theory.   He  specifically  categorizes 
shortcomings  as  being  the  result  of  ignoring  or  overlooking 
one  or  more  of  the  basic  factors  of  production.   First,  the 
model  ignores  the  occurrence  of  price  substitution.   Price 
substitution,  however,  is  a  very  real  phenomenon  since  choice 
among  alternative  inputs  is  very  dependent  on  relative  input 
prices.   Second,  no  accounting  is  made  for  economies  or  dis- 
economies of  scale  for  inputs  which  are  consumed  by  output 
elements.   Third,  certain  factors  of  production,  which  are 
outside  the  model,  are  ignored.   This  criticism  results  from 
a  failure  of  the  model  to  account  for  such  things  as  capital 
stock  depletion  of  natural  resources.   Fourth,  the  model  fails 
to  account  for  joint  production.   This  shortcoming  results 
from  the  establishment  of  a  single  output  measure  for  each  of 
the  model's  input  elements.   An  example  of  this  joint  produc- 
tion is  the  production  of  wool  and  mutton.   During  the  lifetime 
of  a  sheep,  the  hair  from  its  body  is  sheared  and  this  hair 
is  converted  into  wool.   As  the  sheep  matures  and  is  full 
grown,  it  is  slaughtered  and  the  flesh  of  the  animal  is  then 
offered  up  as  mutton.   Thus,  the  sheep  growing  industry  outputs 
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a  joint  product;  wool  from  young  sheep  and  mutton  from  more 
mature  sheep.   Fifth,  and  last,  the  model  fails  to  account 
for  technological  progress  during  the  time-frame  under  con- 
sideration as  the  input-output  coefficients  are  assumed  to 
remain  constant. 

A  second  area  of  criticism  relates  to  weaknesses  that 
may  occur  as  a  result  of  the  manner  in  which  an  actual  model 
is  constructed.   Difficulties  may  arise  as  a  result  of  the 
manner  in  which  industries  are  defined  and  can  be  affected 
by  the  time  when  data  are  measured.   The  result  can  be  that 
the  input-output  coefficients  may  change  as  a  result  of  the 
manner  in  which  the  model  is  aggregated,  or  as  a  result  of  a 
change  in  the  mixture  of  old  and  new  production  procedures . 
To  visualize  the  effect  of  the  manner  of  aggregation,  assume 
that  there  are,  in  fact,  two  constant  coefficients.   If  these 
coefficients  are  added  during  aggregation  of  the  model,  there 
is  no  longer  any  reason  to  assume  that  the  constancy  of  indi- 
vidual coefficients  has  been  preserved.   The  mixing  of  old  and 
new  production  procedures  can  be  clearly  illustrated  by  con- 
sidering an  industrial  firm  that  is  faced  with  a  reduction 
in  demand  for  its  product.   Assume,  for  the  purpose  of  this 
example,  that  an  industry  does,  in  fact,  operate  an  old  and 
a  new  production  process  in  producing  its  output.   When  the 
industry's  sales  decrease,  it  will  undoubtedly  cut  back  on 
production  by  the  old  process.   As  business  improves,  and 
sales  begin  to  increase,  the  firm  will  probably  choose  to 
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expand  by  increasing  its  new  process  and  then  use  the  old 
process  as  capacity  utilization  increases.   The  end  result 
will  be  that  the  relative  proportions  of  old  and  new  produc- 
tion lines  will  have  changed  and  hence,  the  input-output 
coefficients  will  have  changed.   A  final  problem  may  result 
from  the  methodology  used  in  building  the  model.   This  relates 
to  what  Hatanaka  [10]  calls  the  dilemma  of  aggregation.   This 
problem  relates  to  the  two-sided  nature  of  input-output  analy- 
sis, or  the  fact  that  one  must  study  the  industry's  providing 
role  of  furnishing  its  output  as  input  to  an  intermediate 
or  final  output  industry  and  also  the  industry's  consuming 
role  of  using  supplied  inputs  in  the  production  of  outputs. 
The  problem  of  model  validation  can  present  additional 
difficulties.   There  are  essentially  two  approaches  that  may 
be  taken  to  determine  the  ability  of  a  model  to  provide  accurate 
forecasts.   The  difference  in  these  approaches  is  related  only 
to  the  time  period  for  which  the  predictions  are  made.   One 
method  of  prediction  is  to  construct  a  model  now  and  predict 
subsequent  outputs.   The  other  possibility  is  to  use  histori- 
cal data,  available  in  the  economy,  to  build  a  model  and  pre- 
dict outputs  for  which  a  realized  value  already  exists.   "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload"  in  the  11th  Naval  District  was  validated  using  the 
latter  approach. 
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III.   THE  NAVY  PERSONNEL  RESEARCH  AND  DEVELOPMENT 
CENTER'S  REGIONAL  INPUT-OUTPUT  MODEL 

The  Navy  Personnel  Research  and  Development  Center, 
located  at  San  Diego,  California,  using  the  techniques  of 
input-output  modeling, formulated  "A  Regional  Input-Output 
Model  For  Forecasting  Shore-Based  Navy  Workload."   This  model 
was  a  quantitative  tool  to  forecast  workload  at  various  11th 
Naval  District  shore  activities  based  on  a  specific  cleet 
configuration.   The  input-output,  I/O,  model  was  developed  as 
a  means  by  which  I/O  analysis  could  capture  and  accurately 
measure  the  multitude  of  workload  demands  Navy  activities  place 
on  each  other.   It  could  also  help  in  determining  the  feasi- 
bility and  cost  of  developing  a  Navy-wide  I/O  model  and  it 
could  be  used  to  explore  any  mathematical  and/or  computational 
problems  that  might  be  associated  with  the  formulated  model 
[7]  . 

In  developing  a  system  for  allocating  shore-based  manpower 
resources  in  the  11th  Naval  District,  major  emphasis  was  placed 
on  the  design  of  an  I/O  model  to  forecast  the  workload  of 
shore/support  activities  based  on  size  and  configuration  of 
the  fleet.   By  organizing  Navy  shore  activities  into  an  I/O 
matrix,  the  extent  to  which  each  activity  depended  on  every 
other  activity  to  produce  support  could  be  quantified  using 
historical  data.   I/O  analysis  was  seen  as  a  tool  to  determine, 
for  example,  the  impact  on  Long  Beach  Naval  Shipyard's  workload 
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of  introducing  an  additional  destroyer  into  the  overhaul 
schedule.   Also,  it  could  help  estimate  the  increased  work 
necessary  at  the  Naval  Supply  Center,  San  Diego,  California, 
to  service  more  shipyard  requisitions  in  support  of  the  over- 
haul of  the  additional  destroyer.   Thus,  both  direct  and 
indirect  effects  of  the  fleet  were  captured  in  the  model . 

The  regional  I/O  model  was  formulated  as  a  standard  linear 
program  [7:9]  with  the  objective  being  to  minimize  total  activity 
output  while  satisfying  both  final  and  intermediate  demands 
for  that  output.   The  I/O  sectors  represented  11th  Naval  Dis- 
trict shore-based  activities  with  final  demands  represented 

2 

by  the  ships  and  aircraft  of  the  fleet. 

The  regional  I/O  model  formulation  was: 


3 


m 


Minimize    Z   =    T   x. 

iil   ^ 


subject  to: 


(I  -  A)  X  >_  d 
X  >  0 


2 

Shore-based  activities  not  included  as  I/O  sectors  were 

treated  as  final  demands  on  the  system. 

3 
The  I-A  matrix  was  a  square  matrix  and  so  the  solution  to 

the  model  had  to  be  a  point  solution.   Thus,  the  objective 

function  served  no  purpose  except  to  formulate  the  model  as 

a  linear  programming  problem.   Hence,  the  objective  function 

could  be  dropped  since  the  "subject  to"  conditions  determine  the 

unique  solution  to  the  model. 
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where: 

m   =   number  of  I/O  sectors, 

X.   =   total  required  output  level  of  sector  i, 

d.   =   final  demand  on  sector  i, 

A  =  [a. .]  =   amount  of  input  required  by  sector  j  from 

^  sector  i  to  product  one  unit  of  output  from 

sector  j ,  and 

I   =   identity  matrix. 

Final  demand,  d.,  for  a  sector  was  treated  as  a  constant 

1 

on  the  right-hand  side  of  the  linear  programming  formulation. 
Therefore,  calculation  of  final  demand  for  a  sector  by  fleet 
size,  mix,  location,  and  operating  tempo  was  done  outside  the 
LP  model  and  then  summed  up. 

The  formulation  of  the  final  demand  model  was: 


J   K    H 

j 


^^  ^  4i  Ji  h=i  ""^^h  ^^^^  ""i^^^ 


where : 


n.,,    =   number  of  ships  of  type  j,  homeport  k, 
^^^      and  status  h,4 


•jkh 


=   average  time  spent  in  status  h  by  a  ship 
of  type  j  and  homeport  k. 


4 
A  ship  would  be  in  one  of  the  following  statuses:   in 

port,  in  overhaul,  on  local  operations,  or  on  extended  cruise 


24 


r.  .,  ,   =   average  demand  rate  for  sector  i  goods  of 
ships  of  type  j,  homeport  k,  and  status  h, 

J      =   number  of  ship  types, 

K      =   number  of  homeports,  and 

H      =   number  of  statuses  (e.g.,  "in-port"  or 
"over-haul") . 

Nine  major  11th  Naval  District  shore  activities  were 
studied  in  detail  to  develop  an  I/O  model  representing  the 
fleet-support  demand  network.   The  activities  were  selected 
for  their  (1)  wide  range  of  functions,  outputs,  and  data 
problems,  (2)  manpower  intensities,  and  (3)  direct  and  indirect 
linkages  to  the  fleet.   An  I/O  formulation  of  the  fleet-support 
problem  required  the  creation  of  economic  sectors,  each  of 
which  produced  a  single,  unique  output.   While  the  workload 
of  many  shore  activities  could  be  characterized  by  a  single 
output  measure,  other  activities  produced  several  distinct 
outputs.   Thus,  it  was  necessary  to  divide  multioutput  activi- 
ties into  several  sectors,  which  increased  the  number  of  I/O 
sectors  in  the  model  to  twenty.   Demands  placed  on  the  system 
by  units  not  included  in  this  set  of  twenty  sectors  are  all 
treated  as  final  demands  on  the  system.   This  includes  other 
11th  Naval  District  shore  activities.  Navy  activities  in 
other  geographical  areas,  and  the  operating  units,  ships  and 
aircraft,  of  the  fleet.   A  list  of  the  twenty  I/O  sectors  of 
the  model  and  the  units  that  were  used  to  measure  their  total 
workload  are  shown  in  Table  I . 
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Table  I 
I/O  Model  Sectors 


Sector/Activity 


Output  Measure 


1  Naval  Supply  Center 

San  Diego 

2  Long  Beach  Naval  Shipyard 

3  San  Diego-based  Intermediate 

Maintenance  Activities 

4  Naval  Air  Rework  Facility 

North  Island 

5  Public  Works  Center  (PWC) 

San  Diego 

6  PWC — Facilities  Maintenance 


7  PWC — Transportation 

8  PWC — Utilities 

9  Naval  Regional  Medical 

Center,  San  Diego 

10  Naval  Regional  Medical 

Center,  San  Diego 

11  Naval  Station — Military 

Personnel 

12  Naval  Station — Port  Services 

13  NAS  Miramar — Air  Operations 

14  NAS  Miramar — Aircraft 

Intermediate  Maintenance 

15  NAS  Miramar — Supply 

16  NTC — Recruit  Training 

17  NTC — Service  Schools 


Requisitions 


Man-days  of  Ship  Repair 

Man-hours  of  Intermediate 
Maintenance 

Manhours  of  Air  Rework 


Man-hours  of  Direct  Labor 


Man-hours  of  Facilities 
Maintenance 

Gallons  of  Fuel 

Million  BTU's  of  Energy 

Admissions 


Outpatient  Visits 

Military  Personnel  Actions 

Port  Service  Hours 
Air  Operations 
Man-hours  of  AIM 

Requisitions 

Students 

Students 
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Table  I  (Cont'd) 


18  Manpower — Active  Duty  Navy 

19  Manpower — Dependents  of 

Active  Duty  Navy 

20  Manpower — civilian 


Personnel 
Personnel 

Personnel 


A  weighted  average  of  Admissions  and  Outpatient  visits 
was  used  for  the  output  measure. 
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To  develop  an  I/O  coefficient  matrix  for  the  model,  it 
was  necessary  to  study  each  activity  in  the  model.   This  series 
of  investigations  enabled  the  determination  of  each  sector's 
level  of  output  and  the  distribution  of  that  output  among  the 
other  sectors  of  the  model,  as  well  as  the  units  which  make 
up  final  demand  [1; 2; 3; 4 ; 5; 6; 8 ; 12; 14] .   At  each  activity, 
average  demand  rates  were  calculated  over  the  most  recent 
period  of  time  for  which  data  were  available,  generally 
FY75-77. 

In  addition  to  the  flow  of  goods  and  services  between  sec- 
tors of  the  model,  the  analysis  of  activities  that  serve  the 
fleet  directly  required  the  estimation  of  final  demand  rates 
for  the  San  Diego-based  operating  units  of  the  Pacific  Fleet. 
This  involved  examining  the  level  of  demand  placed  on  an 
activity  by  each  of  the  approximately  120  ships  homeported  in 
San  Diego.   These  demand  rates  were  categorized  by  ship  type 
and  status.   Similar  analyses  were  performed  to  determine  the 
level  of  workload  required  to  service  each  type  of  aircraft  at 
the  Naval  Air  Rework  Facility.   Using  these  rates,  the  model 
was  able  to  compute  the  effect  on  final  demand  and  total  out- 
put caused  by  changes  in  the  size  and  deployment  status  of 
operating  units  of  the  fleet. 

Since  workload  at  some  of  the  Navy  activities  could  only 
be  measured  in  terms  of  two  or  more  outputs,  special  estimation 
techniques  had  to  be  used  to  model  the  distribution  of  goods  and 
services  among  model  sectors  that  were  a  part  of  the  same  Navy 
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activity.   It  was  not  feasible  to  attempt  to  trace  each  input 
purchased  by  an  activity  to  the  production  department  from 
which  one  of  its  outputs  came;  hence,  certain  assumptions  were 
made  concerning  the  internal  use  of  inputs  within  these 
activities . 

The  matrix  shown  as  Figure  2  illustrates  the  economic 
interactions  among  sectors  of  the  model.   The  marked  cells 
indicated  the  location  of  nonzero  elements  in  the  I-A  matrix. 
Thus,  a  nonzero  entry  in  element  (i/j)  of  the  matrix  indicated 
that  the  output  of  sector  i  was  used  directly  in  the  production 
of  sector  j's  output.   As  an  example,  element  (1,2)  of  the 
matrix  indicated  an  interaction  between  the  Navy  Supply  Center 
and  Long  Beach  Navay  Shipyard.   This  meant  that  at  least  one 
requisition  was  processed  by  the  Navy  Supply  Center  for  Long 
Beach  Naval  Shipyard  and  that  a  change  in  the  shipyard  workload 
would  in  some  way  impact  the  workload  at  the  supply  center . 

In  an  effort  to  give  the  reader  a  working  knowledge  of 
the  input-output  coefficient  matrix,  the  first  row  of  A  and 
the  first  column  of  A  will  be  derived.   First,  the  input  matrix, 
B,  was  constructed  by  collecting  data  on  inputs  to  each  of  the 
twenty  sectors  of  the  model  from  each  sector  of  the  model .   The 
first  row  of  the  B  matrix  will  now  be  constructed.   Table  II 
is  a  table  of  the  first  row  of  the  input  matrix,  B.   During 
the  base  period  of  FY75-77,  The  Navy  Supply  Center,  NSC, 
generated  6131  requisitions  internally  which  had  to  be  processed 
by  NSC.   Thus,  the  input  to  sector  one,  NSC,  from  sector  one. 
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Figure  2.   11th  Naval  District  Regional  l/O  Model  Impact  Matrix. 
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Table  II 
Row  One  of  the  Input  Matrix,  B 


Sectors   1-  5 

6131 

83387 

98472 

81634 

1090 

Sectors   6-10 

0 

0 

0 

131589 

0 

Sectors  11-15 

0 

6220 

64894 

0 

0 

Sectors  16-20 

9069 

4282 

0 

0 

0 

Read  across  the 

rows. 

NSC,  was  6131  requisitions  and  this  number  was  placed  in  the 
B  matrix  as  b,  , .   Similarly,  Long  Beach  Naval  Shipyard,  sector 
two,  generated  83387  requisitions  which  NSC  processed.   This 
number  was  placed  in  the  B  matrix  as  b,  ^,      Entries  b,  ^,    b,  ., 
and  b,  ^  were  similarly  derived.   All  requisitions  for  the 
Public  Works  Center  were  generated  on  the  administrative  side 
of  the  house  of  the  Public  Works  Center,  hence.  Public  Works 
Center — Facilities  Maintenance,  Public  Works  Center — Trans- 
portation, and  Public  Works  Center — Utilities  had  no  input  of 
requisitions.   Thus,  b,  g,  b^  _,  and  b.  g  were  zero.   The  Naval 
Regional  Medical  Center  at  San  Diego  generated  131589  requisi- 
tions which  were  processed  by  NSC.   Since  requisitions  were 
only  generated  in  the  administrative  side  of  the  house  of  Naval 
Regional  Medical  Center,  all  requisitions  were  put  into  the 
admissions  sector,  sector,  nine,  and  so  the  Naval  Regional 
Medical  Center — Outpatients  sector,  sector  ten,  had  no  input. 
All  requisitions  submitted  to  NSC  from  the  Naval  Station  at 
San  Diego,  6220  requisitions,  were  processed  as  Naval  Station 
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Port  Services  requisitions  and  so  Naval  Station  Military- 
Personnel,  sector  eleven,  had  no  input.   As  a  result  of  this, 
^1  11  ^®^^^^  zero  and  b     became  6220.   NAS  Miramar  sub- 
mitted a  total  of  64894  requisitions  to  NSC  and  all  were  pro- 
cessed as  NAS  Miramar — Air  Operations,  sector  thirteen,  requi- 
sitions.  Thus,  b,  ,-  became  64894.   Since  all  NAS  Miramar 
requisitions  were  processed  as  sector  thirteen  requisitions, 
NAS  Miramar — Aircraft  Intermediate  Maintenace,  sector  fourteen, 
and  NAS  Miramar — Supply,  sector  fifteen,  had  no  inputs  to  the 
matrix.   Consequently,  b,  ,.    and  b,  , -  were  zero.   Naval  Train- 
ing Center — Recruit  Training,  sector  sixteen,  submitted  90  69 
requisitions  to  NSC.   Naval  Training  Center — Service  Schools 
Command,  sector  seventeen,  submitted  4  28  2  requisitions  to  NSC. 
The  last  three  sectors  of  the  model  were  accounting  sectors  to 
track  personnel  only  and  so  no  requisitions  were  submitted  from 
these  sectors.   Thus,  Manpower — Active  Duty  Navy,  sector 
eighteen.  Manpower — Dependents  of  Active  Duty  Navy,  sector 
nineteen,  and  Manpower — Civilian,  sector  twenty,  had  no  inputs 
to  the  B  matrix  for  NSC,  sector  one.   This  completed  the  first 
row  of  the  B  matrix,  the  row  that  represented  the  Naval  Supply 
Center. 

The  first  column  of  the  B  matrix  was  derived  in  a  similar 
manner  to  the  first  row  of  the  matrix  and  Table  III  is  a  table 
of  the  first  column  of  the  input  matrix,  B. 

The  Naval  Supply  Center  input,  b,  ,  remained  unchanged. 
No  man-days  of  ship  repair  from  Long  Beach  Naval  Shipyard, 
sector  two,  were  used  by  NSC,  sector  one,  thus  b^    -^   was  zero. 
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Table  III 
Column  One  of  the  Input  Matrix,  B 


Sectors   1-  5 

6131 

0 

0 

0 

0 

Sectors   6-10 

11734 

69983 

88906 

0 

0 

Sectors  11-15 

0 

0 

0 

0 

0 

Sectors  16-20 

0 

0 

21 

0 

740 

Read  across  the 

rows. 

Similarly,  no  man-hours  of  intermediate  maintenance  from  San 
Diego-based  Intermediate  Maintenance  Activities,  sector  three, 
were  used  by  NSC.   No  man-hours  of  air  rework  from  the  Naval 
Air  Rework  Facility  at  North  Island,  sector  four,  were  used 
by  NSC  and  no  man-hours  of  direct  labor  from  the  Public  Works 
Center,  sector  five,  were  used  by  NSC.   NSC  did  utilize  11734 
man-hours  of  facilities  maintenance  from  the  Public  Works 
Center — Facilities  Maintenance,  sector  six.   The  Public  Works 
Center — Transportation,  sector  seven,  contributed  69983  gallons 
of  fuel  to  the  operation  of  NSC.   In  order  to  heat  buildings 
of  NSC  and  run  auxiliary  equipments.  Public  Works  Center — 
Utilities,  sector  eight,  provided  88906  million  BTU's  of  energy 
to  NSC.   The  only  other  sectors  to  provide  inputs  to  NSC  were 
Manpower — Active  Duty  Navy,  Sector  eighteen,  21  personnel,  and 
Manpower — Civilian,  sector  twenty,  740  personnel.   Since  no 
other  sectors  had  inputs,  bg  ,,  b-Q  ,,  b,,  , ,  b-2  i/  ^^2    i' 

^14,1'  ^15,1'  ^16,1'  ^17,1'  ^^^  ^19,1  ^^^^  ^^^  assigned  the 
value  zero.   This  completed  the  first  column  of  the  B  matrix. 
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The  remainder  of  the  B  matrix  was  filled  in  using  the  same 
techniques  as  just  shown  for  row  one  and  column  one. 

With  the  B  matrix  completed,  the  next  task  was  to  compute 
the  input-output  coefficients  and  put  them  in  the  A  matrix. 
This  was  done  by  dividing  the  inputs  in  the  columns  of  the  B 
matrix  by  the  total  real  output,  during  the  base  period,  of 
the  corresponding  sectors.   Table  IV  is  a  table  of  the  first 
row  of  the  input-output  coefficient  matrix,  A. 

Table  IV 
Row  One  of  the  Input-Output  Coefficient  Matrix,  A 


Sectors   1-  5 

0.00580 

0,10500 

0.03410 

0.01260 

0.00560 

Sectors   6-10 

0.0 

0.0 

0.0 

0.45510 

0.0 

Sectors  11-15 

0.0 

0.62920 

0.59300 

0.0 

0.0 

Sectors  16-20 

1.64000 

0.90800 

0.0 

0.0 

0.0 

Read  across  the 

rows. 

To  illustrate  the  procedure  described  above,  some  of  the 
elements  of  the  input-output  coefficients  for  the  Naval  Supply 
Center,  sector  one,  row  one  of  the  A  matrix  will  now  be  cal- 
culated.  To  compute  a,  ,  element  b   ^  was  divided  by  the  total  real 

J. ,  J.  1  /  j- 

output  for  NSC,  sector  one,  X, .   The  total  real  output  for 
each  of  the  twenty  sectors  may  be  found  in  Table  VI. 


a,  ,   =   b,  ,/X,   =   6131/1057143   =   0.00580 
1,1       1,1^  1 
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The  element  a^  ^  was  computed  by  taking  the  element  b^  ^  and 

^  t  ^  L  ,  2. 

dividing  it  by  the  total  real  output  for  Long  Beach  Naval 
Shipyard,  sector  two,  which  gave 

^1  2  ^  ^1  2^^2  ^  83387/794165  =  0.10500  requisitions/man-day. 

The  element  a^^  ^  was  computed  by  taking  the  element  b,  ^  and 
dividing  it  by  the  total  real  output  for  San  Diego-based 
Intermediate  Maintenance  Activities,  sector  three.   This  gave 

^1  3  ^  ^1  3^^3  ^   98472/2887755  =  0.03410  requisitions/man-hour. 

Since  the  elements  b^^^,  b^^^,  b^^g,  b^^^^^^,  b^^^^,  b^^^^,  b^^^^, 
b,  ,g,  b,  ,Q,  and  b,  20  w^^e  zero  the  elements  a,  g,  a,  _, 

^1,8'  ^1,10'  ^1,11'  ^1,14'  ^1,15'  ^1,18'  ^1,19'  ^^^  ^1,20 
were  assigned  the  value  zero.   This  completed  the  first  row 

of  the  input-output  coefficient  matrix,  A. 

To  calculate  the  values  of  the  input-output  coefficients 

in  the  first  column  of  A,  similar  logic  was  used.   Table  V  is 

a  table  of  the  first  column  of  the  input-output  coefficient 

matrix,  A.   The  first  element,  a,  ,  was  computed  above  and  so 

J-  /  -i- 

it  did  not  have  to  be  recomputed.   The  element  a^  ,  was  com- 
puted by  taking  b^  ,  and  dividing  it  by  the  total  real  output 

b ,  1 

for  NSC,  sector  one,  which  gave 


a^  1  =  b.-  i/X,  =  11734/1057143  =  0.01110  man-hours/requisition. 

D  ,  1       0,1    1 
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Table  V 
Column  One  of  the  Input-Output  Coefficient  Matrix,  A 


Sectors   1-  5 

0.00580 

0.0 

0.0 

0.0 

0.0 

Sectors   6-10 

0.0110 

0.6220 

0.08410 

0.0 

0.0 

Sectors  11-15 

0.0 

0.0 

0.0 

0.0 

0.0 

Sectors  16-20 

0.0 

0.0 

0.00002 

0.0 

0,00070 

Read  across  the 

rows. 

The  element  a2Q  -.  was  computed  by  taking  b^^    ,  and  dividing 
it  by  the  total  real  output  for  NSC,  sector  one,  as 

^20  1  "^  ^20  l^'^l  ^  740/1057143  =  0.00070  personnel/requisition 

The  elements  h^^^,    h^^^,    b^^^,  b3^^,  b^^^^,  b^^^^,  b^^^^, 

^12, r  ^13,1'  ^14,1'  ^15,1'  ^16,1'  ^17,1'  ^^^  ^19,1  ^^^^  ^^^ 
zero.   Thus,  the  corresponding  a. .'s  were  assigned  the  value 

zero. 

When  all  b. .  values  were  divided  by  X.   the  total  real 
output  of  sector  j,  the  complete  input-output  coefficient 
matrix.  A,  was  constructed.   When  the  input-output  coefficient 
matrix  was  subtracted  from  the  identity  matrix,  the  matrix  of 
technological  coefficients,  I-A,  was  generated.   The  techno- 
logical coefficient  matrix  and  the  final  demand  vector  of  the 
model  are  provided  in  Appendix  [A] . 

The  study  which  generated  the  elements  of  the  B  matrix 
also  generated  the  total  real  output  vector  for  the  twenty 

36 


Table  VI 
The  Total  Real  Output  Vector,  X 


Row      Sector 

1  NSC  -  SAN  DIEGO 

2  LBEACH  NAVAL  SHIPYARD 

3  INTER.  MAINT.  ACTIVITY 

4  NAVAIR  REWORK  FACILITY 

5  PUBLIC  WORKS  CENTER 

6  PWC  -  MAINTENANCE 

7  PWC  -  TRANSPORTATION 

8  PWC  -  UTILITIES  ENERGY 

9  NAVREGMEDCEN  -  ADMISS . 

10  NAVREGMEDCEN  -  OUTPAT. 

11  NAVSTA  -  MILITARY  PERS , 

12  NAVSTA  -  PORT  SERVICES 

13  NAS  MIRAMAR  -  AIR  OPS 

14  NAS  MIRAMAR  -  AIMD 

15  NAS  MIRAMAR  -  SUPPLY 

16  NTC  -  RECTRACEN 

17  NTC  -  SERSCCLCCM 

18  MANPOWER  -  ACTIVE  DUTY 

19  MANPOWER  ~  DEPS  OF  ACDU 

20  MANPOWER  -  CIVILIAN 


Output 

1057143.00  REQUISITIONS 

794165.00  MAN-DAYS 

2887755.00  MAN-HOURS 

6478878.00  MAN-HOURS 

194661.00  MAN-HOURS 

1074960.00  MAN-HOURS 

609420.00  GALLONS 

3751626.00  MILLION  BTUS 

289142.00  ADMISSIONS 

498157.00  OUTPATIENTS 

45684.00  PERSONNEL  ACTIONS 

9886,00  PORT  SERVICE  HOURS 

109434,00  AIR  OPERATIONS 

1499711.00  MAN-HOURS 

572336.00  REQUISITIONS 

5530.00  STUDENTS 

4716.00  STUDENTS 

75904.00  PERSONNEL 

17856.00  PERSONNEL 

28209.00  PERSONNEL 
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sectors  of  the  model.   It  is  noteworthy  that  the  row  sums  of 
the  B  matrix  did  not  equal  the  corresponding  elements  of  the 
total  real  output  vector.   This  was  because  the  twenty  sectors 
of  the  model  were  those  of  selected  shore-based  activities  only 
and  the  total  real  output  vector  contained  requirements  for 
the  remaining  shore-based  activities  and  operating  fleet 
units  as  well. 

Using  the  technological  coefficient  matrix  described  above 
and  the  final  demand  vector,  the  solution  to  the  model  was 
found  by  solving 

(I  -A)  X  =  d. 

This  yielded  a  baseline  output  level,  the  output  level  without 
parametric  changes  being  made  to  the  input-output  coefficient 
matrix.   This  was  a  twenty  element  vector  with  the  elements  of 
the  vector  being  the  corresponding  outputs  of  the  sectors  of 
the  model.   This  baseline  output  level  is  presented  as  Table 
VII, 

In  this  section  of  the  thesis,  the  reader  has  been  given 
the  formulation  of  "A  Regional  Input-Output  Model  For  Fore- 
casting Shore-Based  Navy  Workload"  in  the  11  Naval  District 
and  its  component  parts.   The  elements  of  the  input  matrix, 
B,  and  their  origin  have  been  explained.   The  elements  of 
the  input-output  coefficient  matrix.  A,  and  their  computation 
have  been  explained.   The  technological  coefficient  matrix. 
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Table  VII 
Baseline  Output  Levels 


Row 


Sector 


Output  Level 


1  NSC  -  SAN  DIEGO 

2  LB EACH  NAVAL  SHIPYARD 

3  INTER.  MAINT.  ACTIVITIES 

4  NAVAIR  REWORK  FACILITY 

5  PUBLIC  WORKS  CENTER 

6  PWC— MAINTENANCE 

7  PWC — TRANSPORTATION 

8  PWC — UTILITIES  ENERGY 

9  NAVREGMEDCEN — ADMISSIONS 

10  NAVREGMEDCEN— OUTPATIENTS 

11  NAVSTA — MILITARY  PERSONNEL 

ACTIONS 

12  NAVSTA — PORT  SERVICES 

13  NAS  MIRAMAR — AIR  OPERATIONS 

14  NAS  MIRAMAR— AIMD 

15  NAS  MIRAMAR — SUPPLY 

16  NTC — RECRUIT  TRAINING 

CENTER 

17  NTC — SERVICE  SCHOOLS 

COMMAND 

18  MANPOWER— ACTIVE  DUTY  NAVY 

19  MANPOWER — DEPENDENTS  OF 

ACTIVE  DUTY  NAVY 

20  MANPOWER — CIVILIANS 


1572270.00  Requisitions 

794165.00  Man-Days 
2957636.00  Man-Hours 
6478878.00  Man-Hours 
2589904.00  Man-Hours 
1573271.00  Man-Hours 
1503492.00  Gallons  of  Fuel 
5939721.00  Million  BUT's 

310200.00  Admissions 
1726666.00  Outpatients 

45684.00  Personnel  Actions 

9886.00  Port  Service  Hours 
1094  34.00  Air  Operations 
1602596.00  Man-Hours 
620330.00  Requisitions 

5530.00  Students 

4716.00  Students 

106825.00  Personnel 
140704.00  Personnel 

56891.00  Personnel 


I-A,  has  been  presented  and  the  solution  to  the  model,  the 
baseline  output  level,  has  been  presented. 
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IV.   METHODOLOGY  USED  IN  THE  SENSITIVITY  ANALYSIS 

The  Navy  Personnel  Research  and  Development  Center's  "A 
Regional  Input-Output  Model  For  Forecasting  Shore-Based  Navy 
Workload"  in  the  11th  Naval  District  was  composed  of  a  twenty 
by  twenty  technological  coefficient  matrix,  I-A,  and  a  twenty 
element  final  demand  vector,  d.   In  order  to  perform  a  sensi- 
tivity analysis  on  this  model,  it  was  decided  to  parametrically 
change  each  of  the  sixty-three  nonzero  a..'s,  i  =  1,  2,  ..., 
20;  j  =1,  2,  ...,  20,  of  the  input-output  coefficient  matrix, 
A,  by  a  small  percentage  called  delta  and  to  compute  the  output 
of  the  model  with  the  new  technological  coefficient  matrix, 

c  

(I-A*),   and  the  given  final  demand  vector,  d. 

The  elements  of  the  final  demand  vector  were  held  constant 
during  this  process  because  the  interest  was  in  what  would 
happen  to  the  output  of  the  model  if  the  input-output  coeffi- 
cients were  parametrically  changed  for  each  row  of  A  and  a 
solution  was  obtained  to  the  model  with  the  final  demand  vector 
being  held  constant. 

The  decision  was  made  to  change  an  entire  row  rather  than 
to  change  individual  elements  of  the  input-output  coefficient 
matrix.   This  was  because  of  the  manner  in  which  the  a. . 's  were 
computed,  as  discussed  in  Sections  II  and  III  of  this  thesis. 
The  input-output  coefficients  were  computed  from  the  input 


The  asterisk  indicates  a  parametric  change  has  been  effected 
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matrix,  B.   The  elements  of  B  were  taken  from  studies  of  each 
sector  of  the  model,  and  thus,  an  entire  row  of  B  was  the 
result  of  one  study.   Any  one  column  of  B  was  the  partial 
result  of  twenty  such  studies.   It  was  more  likely  that  if 
over-counting  or  under-counting  was  done  in  the  sorting  of 
a  sector's  data  to  get  the  input  for  each  sector,  the  differ- 
ence would  be  more  prevalent  in  the  row  than  in  the  column 
corresponding  to  the  sector  under  study.   Since  the  input 
matrix,  B,  reflected  this  difference  in  its  rows,  it  then 
followed  that  the  a. .'s,  computed  from  the  b. .'s,  would  preserve 
this  property.   From  the  relationship  between  a. .  and  b. ., 
given  in  Section  I,  it  could  be  seen  that  parametrizing  a  row 
of  the  input-output  coefficient  matrix.  A,  by  a  percent  was 
equivalent  to  parametrizing  the  same  row  of  the  input  matrix, 
B,  by  a  percent.   That  is,  since 


a. .   =  b. ./total  real  output  of  industry  j. 


if 


b  .  .   =   a  b  .  . , 


then 


a.  .   =   a  a .  . 
13        ID 


as  long  as  the  total  real  output  of  industry  j  is  not  changed 
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Further,  it  was  found,  on  inspection  and  testing,  that  para- 
metrization  of  rows  produced  consistently  larger  changes  than 
did  parametrization  of  columns.   Thus,  it  was  decided  to  con- 
centrate on  the  effect  that  parametric  changes  to  the  rows 
would  make  and  not  to  concentrate  on  the  effect  that  parametric 
changes  to  the  columns  would  make. 

In  Section  III,  Table  VII,  the  baseline  output  level,  Q, 
of  the  model  was  presented.   Q  was  a  twenty  element  vector 
with  each  element  of  the  vector  representing  one  of  the  twenty 
sectors  of  the  Navy  Personnel  Research  and  Development  Center's 
model.   To  compute  output  levels  of  the  model  with  parametric 
changes  effected  to  the  A  matrix,  each  row  of  the  input-output 
coefficient  matrix  was  parametrically  changed  by  a  percentage, 
one  row  at  a  time,  and  a  new  output  level,  Q  ,  was  obtained. 
One  Q  for  each  row  of  the  input-output  coefficient  matrix 
was  calculated  for  each  delta.   The  elements  of  Q^  represented 
output  levels  for  each  sector  of  the  parametrized  model. 

To  calculate  the  percentage  change  in  output  level  for 
each  of  the  twenty  sectors,  the  baseline  output  level  of  that 
sector  was  subtracted  from  that  sector's  element  in  Q  ,  the 
difference  divided  by  the  baseline  output  level  of  that  sector, 
and  the  result  multiplied  by  100.   The  following  algebraic 
expression  was  utilized  to  compute  the  percentage  changes: 

Q^  -  Q. 
J  ■ — 1  X  100   =  percentage  change 

J 
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in  the  output  of  the  j    sector  of  the  model  when  the  i 
row  of  the  input-output  coefficient  matrix  had  been  para- 
metrized; i  =  1,  2,     ...  ,  20;  j  =  1,  2,    ...  ,  20.   Using  this 
expression,  twenty  percentage  changes  were  computed  for  each 
parametrized  row. 

It  was  necessary  to  decide  which  of  these,  if  any,  were 
significant  enough  to  use  for  later  analysis.   To  do  this,  a 
test-statistic  had  to  be  decided  upon.   It  was  noted  that  most 
percentage  changes  tended  to  be  at  or  near  zero.   It  was  de- 
cided to  use  the  mean  percentage  change  as  the  test-statistic. 
The  mean  percentage  change  for  the  output  of  the  model,  given 
a  parametrized  row,  was  computed.   This  mean  percentage  change 
was  considered  significant  for  the  delta  value  used  to  compute 
the  output  of  theparametrized  model.   Each  percentage  change, 
of  the  twenty  computed  above,  was  then  compared  to  this  test- 
statistic  and  if  it  was  greater  than  the  test-statistic,  it 
was  considered  significant.   The  significant  percentage  changes 
were  identified  and  used  to  further  inspect  the  inter-sectorial 
relationships  in  the  next  section  of  this  thesis.   If  the  per- 
centage change  was  less  than  the  test-statistic,  it  was  con- 
sidered not  significant  and  disregarded. 

As  an  example  of  this  methodology,  significant  percentage 
changes  in  the  output  level  of  the  model  with  row  nine,  the 
row  corresponding  to  Naval  Regional  Medical  Center — Admissions, 
being  changed  by  ten  percent  will  now  be  given.   Table  VIII 
is  the  vector  of  percentage  changes  for  the  model  when  row  nine 
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Table  VIII 

Vector  of  Percentage  Changes  in  Outputs  for  a 
Parametric  Change  in  Row  Nine  of  Ten  Percent 


Sectors 

1-  5 

0.06 

0.0 

0.0 

0.0 

0.03 

Sectors 

6-10 

0.02 

0.02 

0.06 

0.69 

0.11 

Sectors 

11-15 

0.0 

0.0 

0.0 

0.0 

0.0 

Sectors 

16-20 

0.0 

0.0 

0.16 

0.14 

0.13 

Read  across  the  rows;  values  in  percent. 

The  mean  percentage  change  equals  0.07  percent. 


of  the  input-output  coefficient  matrix  was  parametrized  by 

9 
ten  percent.   This  vector  was  written  asQ.,  j=l,  2,    ...  , 

20. 

The  significant  percentage  changes  are  those  greater  than 
the  mean  percentage  change.   Sequentially  going  through  the 
table,  the  significant  percentage  changes  are:   sector  nine, 
0.69  percent;  sector  ten,  0.11  percent;  sector  eighteen,  0.16 
percent;  sector  nineteen,  0.14  percent;  and  sector  twenty, 
0.13  percent. 

Interpretation  of  these  percentage  changes  is  now  easily 
done.   If  the  actual  input-output  coefficients  for  the  Naval 
Regional  Medical  Center — Admissions,  sector  nine,  were  ten 
percent  larger  than  the  estimates,  one  could  expect  the  output 
levels  for  each  of  the  sectors  to  be  the  given  percentages 
greater  than  the  baseline  output  levels.   For  example,  the 
output  for  Manpower — Civilian,  sector  twenty,  would  be  56964 
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or  0.13  percent  greater  than  the  56891  personnel  indicated 
by  the  baseline  output  level. 

The  methodology  employed  in  this  sensitivity  analysis 
was  such  that  a  minus  delta  implied  an  over-estimation  of 
the  input-output  coefficients  by  the  given  percentages  and  a 
positive  delta  implied  an  under-estimation  of  the  input-output 

coefficients.   It  was  reasonable  to  assume  that,  in  the  worst 

7 
case,  delta  should  be  in  the  range  -10  percent  to  +10  percent. 

Thus,  the  computed  percentage  changes  given  in  Section  V  of 
this  thesis  were  for  various  values  of  delta  in  the  range  -25 
percent  to  +25  percent.   This  wide  range  encompassed  the  reason- 
able range  and  allowed  for  a  wider  range  over  which  to  do  the 
analysis. 

Depending  by  how  much  a  particular  row  of  the  input-output 
coefficient  matrix  was  likely  to  be  over-estimated  or  under- 
estimated, lower  or  upper  bounds  on  the  percentage  change  in 
the  output  level  of  the  Navy  Personnel  Research  and  Development 
Center's  model  were  established.   Decreasing  the  coefficients 
of  a  row  by  ten  percent  and  solving  the  model  would  yield  an 
output  level  for  each  of  the  twenty  sectors.   Similarly,  in- 
creasing the  coefficients  of  that  same  row  by  ten  percent  would 
yield  another  output  level  for  each  sector.   The  percentage 
changes  computed  for  each  sector  would  then  be  bounds  on  the 


7 
These  values  were  arrived  at  in  consultation  with  Mr. 

Thomas  Blanco  of  the  Navy  Personnel  Research  and  Development 

Center  who  did  much  of  the  estimation  of  the  input-output 

coefficients. 
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percentage  changes  of  the  respective  sectors.   As  an  example, 
the  coefficients  of  row  nine,  the  row  corresponding  to  the 
Naval  Regional  Medical  Center — Admissions,  were  decreased  by 
ten  percent  and  percentage  changes  in  the  output  levels  of  the 
twenty  sectors  were  calculated.   The  percentage  change  for 
sector  nine  was  -0.69  percent.   When  the  coefficients  were 
increased  by  ten  percent  the  percentage  change  for  sector  nine 
was  0.69  percent.   These  two  percentage  changes  were  lower 
and  upper  bounds  on  the  percentage  change  for  sector  nine  for 
any  values  used  for  parametrization  between  plus  and  minus 
ten  percent.   Decreasing  the  input-output  coefficients  of  row 
nine  by  five  percent  caused  a  percentage  change  of  -0.34 
percent  and  increasing  the  coefficients  of  row  nine  caused  a 
percentage  change  of  0.34  percent.   Both  of  these  percentage 
changes  were  bounded  by  -0.69  percent  and  0.69  percent.   There- 
fore, as  long  as  the  change  in  the  input-output  coefficients 
for  row  nine  was  within  plus  and  minus  ten  percent,  the  resul- 
tant percentage  change  in  the  output  of  sector  nine  would  be 
bound  by  plus  and  minus  0.69  percent.   Further,  the  bounds 
computed  for  the  parametrized  sector  were  also  bounds  on  the 
percentage  change  for  all  other  sectors  of  the  model. 
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V.   RESULTS  OF  THE  PARAMETRIC  SENSITIVITY  ANALYSIS 

In  order  to  conduct  a  parametric  sensitivity  analysis,  a 
computer  algorithm  was  written  following  the  methodology  out- 
lined above  in  Section  IV.   First,  the  algorithm  computed  the 
baseline  output  level ,  the  solution  to  the  Navy  Personnel  Re- 
search and  Development  Center's  Regional  Input-Output  Model 
without  parametrically  changing  the  rows  of  the  input-output 
coefficient  matrix.   Next,  each  row  of  the  input-output  coeffi- 
cient matrix  was  parametrically  changed  by  a  percentage  and  a 
new  output  level  was  calculated.   This  gave  twenty  new  output 
levels,  one  for  each  of  the  twenty  sectors  of  the  model.   Delta 
values  of  -25  percent  to  +25  percent  were  used  in  +5  percent 
increments.   Thus,  a  total  of  ten  values  were  used  for  para- 
metrically changing  each  row  of  the  input-output  coefficient 
matrix.   This  provided  ten  output  levels  for  each  sector  to 
be  compared  with  the  baseline  output  level  of  that  sector. 
The  comparison  was  done  by  computing  the  percentage  change 
from  the  sectors  baseline  output  level  and  identifying  the 
significant  percentage  changes. 

Certain  sectors  of  the  model,  by  the  definition  of  signifi- 
cant percentage  change  given  in  Section  IV,  caused  more  signi- 
ficant percentage  changes  than  others.   Since  the  requirement 
for  manpower  was  a  primary  concern  of  the  Navy  Personnel 
Research  and  Development  Center's  model,  particular  attention 
was  given  to  sectors  which  when  their  input-output  coefficient 
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row  was  changed  caused  significant  percentage  changes  and  to 
all  manpower  sectors.   This  group  included  Naval  Supply  Center-- 
San  Diego,  sector  one;  San  Diego-based  Intermediate  Maintenance 
Activities,  sector  three;  Public  Works  Center,  sector  five; 
PWC — Maintenance,  sector  six;  PWC — Transportation,  sector 
seven;  PWC — Utilities,  sector  eight;  Naval  Regional  Medical 
Center — Admissions,  sector  nine;  Naval  Regional  Medical  Center-- 
Outpatients,  sector  ten;  NAS  Miramar — AIMD,  sector  fourteen; 
NAS  Miramar — Supply,  sector  fifteen;  Manpower — Active  Duty 
Navy,  sector  eighteen;  Manpower — Dependents  of  Active  Duty 
Navy,  sector  nineteen;  and  Manpower — Civilian,  sector  twenty. 
A  summary  of  the  significant  percentage  changes  caused  by 
parametrizing  a  row  of  the  input-output  coefficient  matrix 
corresponding  to  a  particular  sector  is  given  as  Table  IX. 

The  sectors  that  had  no  significant  percentage  changes  were 
ones  whose  row  of  the  input-output  coefficient  matrix  consisted 
of  all  zeros.   This  meant  that  these  sectors  had  no  input 
from  the  other  sectors  of  the  model.   Thus,  their  rows  of 
the  input  matrix  were  all  zero.   Since  the  input  matrix  rows 
and  the  input-output  coefficient  matrix  rows  were  all  zero, 
on  parametrization  the  coefficients  remained  zero.   Solution 
of  the  parametrized  model  then  gave  the  same  output  levels  as 
the  baseline  output  level  and  thus  no  significant  percentage 
changes  were  computed.   The  outputs  of  these  sectors  were 
always  constant. 

Line  numbering  on  the  figures  of  this  section  was  done  to 
indicate  the  sector  the  line  corresponded  to.   The  lines  were 
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Table  IX 

Summary  of  Significant  Percentage  Changes  Caused  by 
Parametrizing  a  Particular  Sector 


Sector  Parametrized 


Sectors  with  Significant 
Percentage  Changes 


1  NAVAL  SUPPLY  CENTER 

2  LBEACH  NAVAL  SHIPYARD 

3  INTER.  MAINT.  ACTIVITIES 

4  NAVAIR  REWORK  FACILITY 

5  PUBLIC  WORKS  CENTER 

6  PWC— MAINTENANCE 

7  PWC — TRANSPORTATION 

8  PWC — UTILITIES  ENERGY 

9  NAVREGMEDCEN — ADMISSIONS 

10  NAVREGMEDCEN — OUTPATIENTS 

11  NAVSTA — MILITARY  PERSONNEL 

12  NAVSTA — PORT  SERVICES 

13  NAS  MIRAMAR — AIR  OPS 

14  NAS  MIRAMAR — AIMD 

15  NAS  MIRAMAR — SUPPLY 

16  NTC— RECTRACEN 

17  NTC — SERSCOLCOM 

18  MANPOWER— ACTIVE  DUTY  NAVY 

19  MANPOWER — DEPENDENTS  OF  ACDU 

20  MANPOWER — CIVILIAN 


1,  7 
None 
1/  3 
None 
5,  7 
5,  6,  7 
5,  7 
5,  8 

9,  10,  18,  19,  20 
10 

None 

None 

None 

14,  15 

15 

None 

None 

10,  18,  19 
10,  19 

20 


provided  so  the  reader  could  have  a  visual  representation  of 
the  percentage  change  computed  for  a  known  value  by  which  the 
input-output  coefficient  matrix  row  for  the  sector  in  the  figure 
heading  was  changed.   Further,  to  visualize  the  bounds,  the 
reader  may  select  a  percentage  change  for  a  sector  of  concern. 
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The  bound  will  be  the  ordinate  of  the  least-squares  line  at 
the  point  on  the  line  corresponding  to  the  chosen  percentage 
change . 

The  least-squares  lines  on  some  figures  were  very  flat, 
indicating  small  slopes.   This  meant  that  the  percentage 
changes  resulting  from  changes  in  the  values  of  the  rows  of 
the  input-output  coefficient  matrix  were  very  small.   Some 
of  these  "flat"  lines  may  seem  uninteresting  to  the  reader 
but  they  were  included  for  the  many  managers  who  may  read 
this  thesis,   A  quick  look  at  a  graph  is  sometimes  easier  to 
understand  than  to  have  to  compute  values  and  do  the  graphing 

for  one's  self. 

2 

The  coefficient  of  determination,  r  ,  was  computed  to 

establish  how  well  the  delta  values  and  their  associated 

percentage  changes  for  a  significant  sector  fit  the  linear 

2  2 

regression.   AT  r   equal  zero,  you  have  no  fit  and  at  r   equal 

one,  you  have  a  perfect  fit.   In  the  case  of  all  sectors  which 

2 
showed  a  significant  percentage  change,  r  was  approximately 

2 
equal  to  one.   The  r   values  were  m  all  cases  greater  than 

0.9995.   Thus,  the  percentage  changes  computed  for  the  given 
delta  values  were  very  close  to  being  linear.   Since  the  model 
was  linear  and  the  parametric  changes  formed  linear  transfor- 
mations of  the  input-output  coefficient  rows,  the  resultant 
percentage  changes,  computed  from  the  sectorial  outputs,  would 
themselves  be  expected  to  be  linear.   It  was  felt  that  the 
reason  the  regression  line  did  not  pass  exactly  through  the 
origin  was  because  of  round-off  by  the  computer  in  calculating 
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the  inverse  of  the  twenty  technological  matrices,  I-A* .   Some 
of  the  percentage  changes  were  not  considered  to  be  significant 
and  so  the  regression  lines  for  their  sectors  were  not  graphed. 
Nonetheless,  the  reader  might  be  interested  in  the  magnitude 
of  the  nonsignificant  percentage  changes  so  all  percentage 
changes  resulting  from  parametrization  of  a  model  sector  are 
given  in  Appendix  [B] .   The  reader  may  check  this  appendix 
for  any  sector  of  concern. 

A.   NAVAL  SUPPLY  CENTER 

In  terms  of  significant  percentage  changes,  for  all  delta 
values  used  in  the  sensitivity  analysis,  parametric  changes 
to  the  input-output  coefficients  of  sector  one.  Naval  Supply 
Center — San  Diego,  produced  significant  changes  in  the  output 
levels  of  Navl  Supply  Center,  sector  one,  and  PWC — Transporta- 
tion, sector  seven.   This  could  be  reasoned  out  since  in  order 
for  the  Naval  Supply  Center  to  serve  the  other  sectors,  supply 
them  with  goods,  the  Naval  Supply  Center  must  have  transporta- 
tion to  transport  goods  from  the  supply  depot  to  the  requisi- 
tioning activity.   Thus,  the  significant  change  in  the  output 
level  of  PWC — Transportation,  given  the  change  in  the  Naval 
Supply  Center's  output  was  expected.   Other  sectors  showed 
change  as  well,  but  based  on  the  definition  of  "significant 
percentage  change,"  the  sectors  of  concern  were  the  Naval 
Supply  Center  and  PWC — Transportation.   The  significant  percen- 
tage changes  that  occurred  as  a  result  of  parametrizing  row 
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one  of  the  input-output  coefficient  matrix  are  graphed  as 
Figure  3. 

B.  INTERMEDIATE  MAINTENANCE  ACTIVITIES 

Parametric  changes  to  the  nonzero  elements  of  the  input- 
output  coefficients  in  row  three,  the  row  corresponding  to 
San  Diego-based  Intermediate  Maintenance  Activities,  produced 
significant  changes  in  the  output  levels  of  Naval  Supply 
Center — San  Diego,  sector  one,  and  San  Diego-based  Intermediate 
Maintenance  Activities.   Knowing  that  the  Intermediate  Main- 
tenance Activities  were  heavily  dependent  on  the  Naval  Supply 
Center  to  supply  repair  parts  for  maintenance,  this  linkage 
was  to  be  expected.   The  significant  percentage  changes  that 
resulted  from  parametrizing  row  three  are  graphed  as  Figure  4. 

C.  PUBLIC  WORKS  CENTER 

Changing  the  nonzero  elements  of  row  five  of  the  input- 
output  coefficient  matrix,  the  row  representing  the  Public 
Works  Center,  sector  five,  produced  significant  changes  in 
the  output  levels  of  Public  Works  Center,  sector  five,  and 
PWC — Transportation,  sector  seven.   Most  of  the  output  for  the 
Public  Works  Center  required  transportation  of  some  sort  to 
expend  the  man-hours  so  the  change  in  the  output  level  of 
PWC — Transportation  was  to  be  expected.   One  would  also  expect 
changes  in  the  output  levels  of  the  sectors  which  supported 
the  Public  Works  Center,  i.e..  Naval  Supply  Center — San  Diego, 
PWC — Utilities,  Manpower — Civilian,  and  as  expected,  changes 
occurred.   However,  by  the  definition  of  "significant,"  these 
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changes  were  not  significant.   It  was  also  noted  that  the 
percentage  change  in  output  of  the  Public  Works  Center  was 
approximately  proportional  to  the  amount  of  change  made  to 
the  input-output  coefficients  of  that  sector.   Figure  5  is 
a  graph  of  the  significant  percentage  changes  resulting  from 
parametrizing  row  five. 

D  .   PWC — MAINTENANCE 

Parametrically  changing  the  nonzero  elements  of  row  six 
of  the  input-output  coefficient  matrix,  the  row  corresponding 
to  PWC — Maintenance,  sector  six,  produced  significant  changes 
in  the  output  levels  of  Public  Works  Center,  sector  five, 
PWC--Maintenance,  sector  six,  and  PWC — Transportation,  sector 
seven.   Since  the  model  was  originally  formulated  with  all  PWC 
sectors  combined  into  the  one  sector.  Public  Works  Center, 
the  change  in  the  output  level  of  the  Public  Works  Center  was 
expected.   Since  maintenance  was  required  to  keep  the  fiscal 
plants  of  the  Public  Works  Center  up,  the  change  in  the  output 
level  of  sector  six,  PWC — Maintenance,  was  expected.   The 
change  in  PWC — Transportation  was  expected  also  since  trans- 
portation was  required  to  transport  personnel  and  equipment 
to  maintenance  tasks  in  the  11th  Naval  District.   Figure  6  is 
the  graph  of  significant  percentage  changes  caused  by  para- 
metrizing row  six. 

E.   PWC — TRANSPORTATION 

Changing  the  elements  of  row  seven,  the  row  corresponding 
to  sector  seven,  PWC — Transportation,  evidenced  significant 
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changes  in  the  output  levels  of  Public  Works  Center,  sector 
five,  and  PWC — Transportation,  sector  seven.   Inasmuch  as 
the  relationship  between  all  sectors  of  the  Public  Works 
Center,  as  explained  above,  existed,  these  changes  were 
expected.   The  graph  of  the  significant  percentage  changes 
resulting  from  parametrizing  row  seven  is  presented  as  Figure  7. 

F.  PWC — UTILITIES 

Changing  row  eight,  the  row  corresponding  to  PWC — Utilities, 
sector  eight,  caused  significant  changes  to  Public  Works 
Center,  sector  five,  and  PWC — Utilities,  sector  eight.   As 
explained  above,  the  relationship  of  all  Public  Works  Center 
sectors  would  lead  one  to  expect  these  changes.   The  graph  of 
these  significant  percentage  changes  is  given  as  Figure  8. 

G.  NAVREGMEDCEN — ADMISSIONS 

The  sector  whose  input-output  coefficient  row,  when 
changed,  caused  the  greatest  number  of  significant  percentage 
changes  was  the  Naval  Regional  Medical  Center — Admissions, 
sector  nine.   Changing  the  nonzero  elements  of  row  nine  pro- 
duced significant  changes  in  the  output  levels  of  Naval  Regional 
Medical  Center — Admissions,  Naval  Regional  Medical  Center — 
Outpatients,  Manpower — Active  Duty  Navy,  Manpower — Dependents 
of  Active  Duty  Navy,  and  Manpower — Civilian.   The  change  in 
the  output  level  of  Naval  Regional  Medical  Center — Admissions 
was  to  be  expected  since  the  input-output  coefficients  in  that 
row  were  changed.   A  change  in  the  number  of  personnel  admitted 
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to  the  Naval  Regional  Medical  Center  would  generally  cause 
a  change  in  the  number  of  outpatient  visits  and  so  a  change 
in  outpatient  visits  was  to  be  expected.   Given  that  active 
duty  personnel  and  dependents  of  active  duty  personnel  com- 
prise a  large  part  of  the  number  of  patients  who  are  seen  at 
the  hospital,  change  in  the  output  of  the  sectors  corresponding 
to  these  two  groups  of  people,  sectors  eighteen  and  nineteen 
respectively,  was  expected.   Civilians  make  up  a  large  part 
of  the  staff  at  the  Naval  Regional  Medical  Center  and  so 
change  in  the  workload  at  the  hospital  would  precipitate  a 
change  in  the  civilian  workforce  at  the  hospital.   These  signi- 
ficant percentage  changes  are  given  as  the  graph  of  Figure  9. 

H.   NAVREGMEDCEN — OUTPATIENTS 

Changing  row  ten  of  the  input-output  coefficient  matrix, 
the  row  which  corresponded  to  Naval  Regional  Medical  Center — 
Outpatients,  sector  ten,  showed  only  a  change  in  the  output 
of  that  sector.   The  Naval  Regional  Medical  Center  was  originally 
only  one  sector  in  the  model  and  so  it  was  assumed  that  change 
in  the  output  level  of  the  related  sector  was  absorbed  in  the 
Naval  Regional  Medical  Center — Admissions  sector,  sector  nine. 
Figure  10  is  the  graph  of  significant  percentage  changes. 

I.   NAS  MIRAMAR — AIMD 

Changing  row  fourteen,  the  row  corresponding  to  NAS 
Miramar — AIMD,  sector  fourteen,  produced  significant  changes 
in  the  outputs  of  NAS  Miramar — AIMD,  sector  fourteen,  and  NAS 
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Miramar--Supply ,  sector  fifteen.   This  was  to  be  expected 
since  the  Aircraft  Intermediate  Maintenance  Depot  is  heavily 
dependent  on  NAS  Miramar — Supply  for  repair  parts  to  effect 
repairs  to  aircraft.   See  Figure  11  for  a  graph  of  the  signi- 
ficant percentage  changes. 

J.   NAS  MIRAMAR — SUPPLY 

Changing  row  fifteen,  the  row  corresponding  to  NAS  Miramar — 
Supply,  sector  fifteen,  produced  significant  change  in  the 
output  level  of  NAS  Miramar — Supply,  sector  fifteen,  only. 
This  was  to  be  expected  since  all  requisitions  from  NAS 
Miramar  were  processed  by  NSC  as  NAS  Miramar--Air  Operations, 
requisitions.   The  percentage  changes  ranged  from  -1.94  percent 
for  a  delta  value  of  -25  percent  to  +1.94  percent  for  a  delta 
value  of  +25  percent.   Since  the  sector  under  study  was  the 
only  sector  to  show  significant  percentage  change,  no  graph 
is  provided. 

K.   MANPOWER — ACTIVE  DUTY  NAVY 

Parametrically  changing  the  nonzero  elements  of  row  eighteen, 
the  row  corresponding  to  Manpower — Active  Duty  Navy,  sector 
eighteen,  produced  significant  percentage  changes  in  the  output 
levels  of  Naval  Regional  Medical  Center--Outpatients ,  sector 
ten.  Manpower — Active  Duty  Navy,  sector  eighteen,  and  Man- 
power— Dependents  of  Active  Duty  Navy,  sector  nineteen.   Changing 
the  number  of  active  duty  personnel  in  the  11th  Naval  Dis- 
trict should  cause  a  change  in  the  number  of  dependents  of 
active  duty  personnel.   With  a  change  in  the  number  of  active 
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duty  personnel  and  a  change  in  the  number  of  dependents  of 
active  duty  personnel  in  the  11th  Naval  District  one  would 
expect  a  change  to  occur  in  the  number  of  outpatient  visits 
at  the  Naval  REgional  Medical  Center.   Significant  percentage 
changes  are  graphed  as  Figure  12. 

L.   MANPOWER — DEPENDENTS  OF  ACTIVE  DUTY 

Changing  row  nineteen,  the  row  corresponding  to  Manpower — 
Dependents  of  Active  Duty  Navy,  sector  nineteen,  produced 
significant  changes  in  the  output  levels  of  Naval  Regional 
Medical  Center — Outpatients,  sector  ten,  and  Manpower- -Dependents 
of  Active  Duty  Navy,  sector  nineteen.   This  was  to  be  expected 
since  dependents  of  active  duty  personnel  comprise  a  large 
part  of  the  Naval  Regional  Medical  Center's  outpatient  service. 
See  Figure  13  for  the  graph  of  significant  percentage  changes. 

M.   MANPOWER — CIVILIAN 

Changing  row  twenty  of  the  input-output  coefficient  matrix, 
the  row  corresponding  to  Manpower — Civilian,  sector  twenty, 
produced  a  significant  change  in  the  output  level  of  Manpower 
Civilian,  sector  twenty.   One  would  expect  changes  to  occur 
in  the  output  levels  of  all  sectors  where  civilian  manpower 
was  utilized  but  the  model  only  showed  a  change  of  any  kind 
in  the  output  for  Manpower — Civilian.   Figure  14  is  the  graph 
of  this  significant  percentage  change. 

As  discussed  above,  upper  and  lower  bounds  may  be  taken 
from  the  graphs.   Choose  the  plus  and  minus  percentages  a  row 
is  most  likely  to  be  off  by.   Then  construct  the  necessary 
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Figure   12-     MANPOWER— Active  Duty  Navy,   Sector  18. 
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Figure  13-   MANPOWER— Dependents  of  Active  Duty,  Sector  19* 
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perpendiculars  from  the  chosen  values  to  the  least-squares 
line  of  the  sector  of  interest.   From  the  points  of  inter- 
section on  the  least-squares  line  to  the  ordinate  axis  of 
the  graph  construct  other  perpendiculars.   Where  the  normal 
lines  intersect  the  ordinate  axis,  bounds  may  be  read. 
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VI.   SUMMARY 

A.   CONCLUSIONS 

This  section  presents  the  conclusions  resulting  from  the 
parametric  sensitivity  analysis. 

The  most  sensitive  sector  to  parametric  change  of  its 
input-output  coefficients  was  the  Public  Works  Center.   When 
its  input-output  coefficient  row,  row  five,  was  parametrized 
and  percentage  changes  were  computed  for  the  various  values 
of  delta,  this  sector  had  a  least-squares  line  with  a  slope 
of  0.98.   This  implied  that  if  row  five,  the  row  of  the  A 
matrix  corresponding  to  the  Public  Works  Center  was  parametrized, 
the  percentage  change  in  the  output  of  that  sector  would  be 
approximately  equal  to  the  percentage  used  for  parametriza- 
tion.   Coupled  with  this  was  the  percentage  change  of  PWC — 
Transportation.   The  slope  of  the  least-squares  line  for  PWC — 
Transportation,  sector  seven,  was  0.23.   This  implied  that  if 
the  row  of  the  input-output  coefficient  matrix  representing 
the  Public  Works  center  were  changed  by  a  small  amount  then 
the  output  level  of  PWC — Transportation  would  be  approximately 
that  of  the  baseline  output  level  of  PWC — Transportation  plus 
this  baseline  output  level  times  0.23  times  delta.   This  inter- 
sectorial  relationship  indicated  that  the  managers  of  these 
two  sectors  would  want  to  be  very  mindful  of  changes  in  the 
output  level  of  the  Public  Works  Center  as  a  change  in  its 
output  level  would  precipitate  a  significant  percentage 
change  in  the  output  level  of  PWC — Transportation. 
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The  second  most  sensitive  sector  was  Manpower — Dependents 
of  Active  Duty,  sector  nineteen.   This  sector  had  a  least- 
squares  line  with  a  slope  of  0.88.   Thus,  changing  its  input- 
output  coefficient  row  would  cause  the  output  level  of  the 
model,  with  the  parametrized  row  to  be  approximately  that  of 
the  baseline  output  level  for  sector  nineteen  plus  the  baseline 
output  level  times  0,88  times  delta.   The  intersectorial 
relationship  between  this  sector  and  Naval  Regional  Medical 
Center — Outpatients,  sector  ten,  was  such  that  changing  row 
nineteen  of  the  A  matrix  produced  significant  percentage 
changes  in  the  output  level  of  Naval  Regional  Medical  Center — 
Outpatients.   Further,  changing  the  coefficients  representing 
Manpower- -Dependents  of  Active  Duty  produced  a  Naval  Regional 
Medical  Center — Outpatients  least-squares  line  with  a  slope 
of  0.27.   Hence,  changing  row  nineteen  would  cause  the  output 
of  sector  ten  to  be  approximately  that  of  the  baseline  output 
level  for  sector  ten  plus  the  baseline  output  level  times 
0.27  times  delta.   Thus,  a  change  in  the  number  of  dependents 
would  precipitate  a  change  in  the  number  of  outpatients  who 
visited  the  Naval  Regional  Medical  Center. 

The  third  most  sensitive  sector  was  Naval  Regional  Medi- 
cal Center — Outpatients,  sector  ten.   For  this  sector,  a 
least-squares  line  slope  of  0.71  was  computed.   This  implied 
that  if  the  input-output  coefficients  of  row  ten  were  changed, 
the  output  of  the  model  with  the  changed  coefficients  would 
be  approximately  equal  to  the  baseline  output  level  for  sector 
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ten  plus  the  baseline  output  level  for  sector  ten  times  0.71 
times  delta. 

The  fourth  most  sensitive  sector  was  PWC — Transportation, 
sector  seven.   Changing  the  input-output  coefficients  of  row 
seven  yielded  a  line  with  a  slope  of  0.62  when  the  output 
level  of  the  model,  with  the  changed  coefficients,  was  com- 
pared to  the  baseline  output  level  of  sector  seven.   Thus, 
changing  the  coefficients  of  row  seven  would  mean  that  the 
output  of  sector  seven  would  be  approximately  equal  to  the 
baseline  output  level  of  sector  seven  plus  the  baseline  out- 
put level  of  sector  seven  times  0.62  times  delta.   Coupled 
with  this  was  a  significant  percentage  change  in  the  output 
level  of  Public  Works  Center,  sector  five.   Making  the  afore- 
mentioned changes  to  row  seven  would  produce  a  least-squares 
line  for  the  Public  Works  Center  with  a  slope  of  0.19.   There- 
fore, by  changing  row  seven  the  output  of  the  Public  Works 
Center  would  be  approximately  equal  to  the  baseline  output 
level  of  the  Public  Works  Center  plus  that  baseline  output 
level  times  0.19  times  delta.   In  terms  of  number  of  man-hours, 
this  number  was  not  too  great,  however,  by  the  definition  of 
significant  percentage  change,  it  was  significant. 

The  fifth  most  sensitive  sector  to  parametrization  was 
Manpower — Civilian,  sector  twenty.   Changing  the  coefficients 
in  row  twenty  produced  a  line  with  a  slope  of  0.50  when  the 
output  level  of  the  model,  with  the  parametrized  coefficients, 
for  sector  twenty  was  compared  to  the  baseline  output  level 
for  sector  twenty.   Thus,  changing  row  twenty  by  an  amount 
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delta  would  produce  an  output  level  approximately  equal  to 
the  baseline  output  level  for  Manpower — Civilian,  sector 
twenty,  plus  this  baseline  output  level  times  0.50  times 
delta. 

There  were  other  sectors  one  could  concentrate  on,  how- 
ever, all  others  had  least-squares  line  slopes  of  less  than 
0.50.   The  reader  might  check  the  others  by  going  through  the 
figures  of  Section  V  and  comparing  the  slopes  of  the  least- 
squares  lines.   The  measure  of  sensitivity  used  to  determine 
the  most  sensitive  sectors  was  the  slope  of  the  least-squares 
line  through  the  points  generated  by  computing  percentage 
changes  at  the  given  values  of  delta. 

The  relationship  of  a  change  in  a  row  of  the  input-output 
coefficient  matrix  to  a  change  in  the  original  input  matrix 
was  a  simple  one.   The  input-output  coefficients,  a..'s,  were 
computed  from  the  input  matrix,  B,  and  the  total  real  output 
vector,  X.   The  algebraic  expression  for  this  transformation 
was 


a.  .   =  b.  ./x. 
i:      ID  3 


and  since  elements  of  the  total  real  output  vector  were  held 
constant,  it  was  simple  to  compute  the  inverse  of  this  trans- 
formation given  a  change  to  a. .   To  show  the  inverse  transfor- 
mation, assume  row  i  of  the  input-output  coefficient  matrix 
was  changed  by  an  amount  k.   Then 
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(l+k)a.    =   b.  */x. 


and  solving  this  expression  for  b.   yielded 


b.  *   =   (1  +k)a.  X.  . 
1*  !•  3 


Now,  let's  assume  the  elements  of  row  one  of  the  A  matrix 

* 

were  decreased  by  ten  percent  and  compute  new  values  of  b,  , 

* 

where  b,  ,  would  be  the  value  of  b,  ,  had  the  original  b, 
1,1  1,1  ='1* 

been  such  that  the  values  of  a, ^*  were  computed  from  them. 
Then , 


(l-O.lO)a^^^   =  ^l,l/^l 


which  on  substitution  would  be 


(1  -  0.10)  (0.00580)   =  b,  ,/1057143  requisitions 


or 


(0.90) (0.00580)   =  b,  ^71057143  requisitions 
Solving  this  would  yield 


b,  ,   =   5518  requisitions 
1,1 
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since  5518  is  90  percent  of  6131,  the  original  b,  , ,  this 
shows  that  a  change  of  k  to  the  a. ^  implied  a  change  of  k 
to  the  corresponding  b.  .   Simply  stated,  changing  a.   by 
some  amount  yielded  a  corresponding  change  of  the  same  amoiint 
to  the  respective  b.  . 

The  results  of  this  analysis  may  be  used  to  evaluate  a 
sector's  output  based  on  a  historical  change  in  the  sector's 
performance.   Suppose  a  sector's  actual  inputs  from  the  other 
sectors  were  increased  by  some  percentage.   The  model's  output, 
when  the  model  was  changed  to  reflect  the  new  inputs,  could  be 
predicted  by  finding  the  amount  of  the  change  on  the  appropriate 
figure  or  in  Appendix  [B]  and  reading  the  percentage  change 
from  the  baseline  output  level. 

Suppose  managers  of  the  Naval  Supply  Center  and  the  Naval 
Regional  Medical  Center  were  to  ask  themselves  "what  if  the 
model's  input-output  coefficients  were  off  by  plus  and/or  minus 
ten  percent?   What  range  of  requisitions  or  admissions  might 
the  respective  sectors  encounter."   Thus,  the  sectors  managers 
would  be  interested  in  total  variability  of  forecasts  with 
the  given  model  input-output  coefficients. 

If  the  inputs  to  the  Naval  Supply  Center  were  actually  ten 
percent  greater  in  magnitude  than  the  inputs  in  the  original 
model,  the  percentage  change  in  the  output  of  the  Naval  Supply 
Center  would  be  3.30  percent.   This  would  imply  that  instead 
of  an  output  of  1572270  requisitions  the  Naval  Supply  Center 
would  have  an  output  of  1624129  requisitions.   Similarly,  for 
the  same  change  in  the  input  so  the  Naval  Supply  Center  the 
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output  of  PWC — Transportation  would  be  0.25  percent  greater 
than  the  baseline  output.   The  output  of  PWC — Transportation 
would  increase  from  1503492  gallons  of  fuel  to  1507196  gallons 
of  fuel.   Thus,  in  order  for  the  Naval  Supply  Center  to  output 
an  additional  51859  requisitions,  the  Naval  Supply  Center 
would  precipitate  the  use  of  an  additional  3704  gallons  of 
fuel  from  PWC — Transportation.   Further  the  Naval  Supply  Cen- 
ter's change  in  output  would  cause  a  change  in  the  output  of 
Manpower--Civilian  of  0.07  percent  to  56929  manhours  of  an 
additional  38  manhours  of  civilian  labor. 

If  the  inputs  to  the  Naval  Supply  Center  were  ten  percent 
less  than  the  inputs  in  the  original  model,  the  percentage 
change  in  the  output  of  the  Naval  Supply  Center  would  be 
-3.29  percent.   This  would  imply  an  output  of  1520477  requi- 
sitions.  The  output  of  PWC — Transportation  would  be  1499794 
gallons  of  fuel  instead  of  1503492  gallons  of  fuel,  a  change 
of  -0.25  percent.   The  output  of  Manpower — Civilian  would  be 
56853  manhours  instead  of  56890  manhours,  a  change  of  -0.07 
percent.   The  above  figures  would  be  ranges  on  the  output  levels 
of  significant  sectors,  given  that  the  coefficients  for  the 
Naval  Supply  Center — San  Diego  were  within  plus  and  minus  ten 
percent.   These  ranges  would  be:   NSC — San  Diego  1520477 
to  1624129  requisitions;  PWC — Transportation  1499794  to 
1507196  gallons  of  fuel;  and  Manpower — Civilian  56853  to 
56929  manhours.   These  ranges  would  reflect  both  the  direct 
and  indirect  inputs  required  for  the  Naval  Supply  Center  to 
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operate  with  input-output  coefficients  within  plus  and  minus 
ten  percent  of  the  original  coefficients  of  the  model. 

The  administrator  of  the  Naval  Regional  Medical  Center 
might  wonder  what  services  would  be  needed  if  the  number  of 
shore-based  military  personnel  in  the  11th  Naval  District  was 
increased  or  decreased  by  ten  percent.   The  administrator  would 
only  have  to  go  to  Appendix  [B]  and  look  under  the  ten  percent 
column  of  Manpower — Active  Duty  and  read  off  the  percentage 
changes  for  the  two  sectors  of  the  Naval  Regional  Medical  Center 
and  see  that  there  would  be  0.19  percent  or  581  more  admissions 
for  a  total  of  310782  admissions  and  1.98  percent  or  34207 
more  outpatient  visits  for  a  total  of  1760873  outpatient 
visits.   The  minus  ten  percent  column  would  reveal  a  -0.19 
percent  change  or  an  output  of  309621  admissions  and  a  -1.98 
percent  change  or  an  output  of  1692562  outpatient  visits. 
ThuS/  the  range  of  outputs  would  be:   Naval  Regional  Medical 
Center — Admissions  309621  to  310782  admissions;  and  Naval 
Regional  Medical  Center — Outpatients  1692562  to  1760873  out- 
patient visits. 

Managers  of  the  various  sectors  of  this  model  might  use 
these  results  to  predict  future  requirements  for  their  sectors. 
It  was  pointed  out  above  that  predictions  could  be  made  as  to 
increased  requirements  based  on  increased  outputs  for  one  sec- 
tor.  Similarly,  decreased  requirements  and/or  resources  would 
result  from  a  decrease  in  the  input  to  a  sector  which  would 
necessitate  a  decrease  in  the  output  of  that  sector. 
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B.   RECOMMENDATIONS 

Two  specific  areas  for  further  study  are  recoimnended. 

First,  the  model's  input-output  coefficient  matrix  needs 
to  be  reevaluated  to  see  if  all  the  zero  entries  are  still 
zero.   For  instance,  when  model  sectors  of  the  original  model 
formulation  were  split  to  form  other  sectors  were  the  inputs 
accurately  and  correctly  recorded?  One  would  think  that  per- 
haps there  were  requisitions  submitted  to  the  Naval  Supply 
Center  in  support  of  outpatient  visits  to  the  Naval  Regional 
Medical  Center,  but  none  existed  in  the  model  as  formulated. 
Perhaps  it  was  not  possible  to  separate  the  supplies  that  went 
to  the  Naval  Regional  Medical  Center  because  of  the  manner  in 
which  accounting  was  done.   Nonetheless,  this  is  an  area  for 
further  study,  if  possible,  or  perhaps  a  means  of  tracking 
supplies  used  by  both  sectors  of  the  Naval  Regional  Medical 
Center  might  be  initiated  if  it  has  not  already  been  done. 

The  second  area  for  possible  further  study  is  all  sectors 
of  the  Public  Works  Center.   The  Public  Works  Center  was  the 
most  sensitive  sector  of  the  model  and  divisions  of  the  PL±)lic 
Works  Center  were  all  related  so  changes  in  the  coefficient 
matrix  here  could  be  crucial  to  the  model's  output.   Further 
work  in  this  area  could  prove  to  be  well  worth  the  effort  in 
terms  of  accurate  forecasts  by  the  model. 
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The  Final  Demand  Vector »  d 


Row       Sector 

1  NSC  -  S/iN  CIEGC 

2  LBEACh  hAV4L  SHIPVARC 

3  INTER,  VAIM.  ACTIVITY 
^  NAVAIF  PEUCRK  FACILITY 

5  PUBLIC  WCPHS  CENTER 

6  PWC  -  MAINTENANCE 

7  PVC  -  TRANSPQRTATICN 

8  FWC  -  L7UITIES  ENERGY 

9  NAVREGMEDCEN  -  AOMISS, 

10  NAVREGKEOCEN  -  CUTPAT, 

11  N)!\STA  -  MILITARY  FERS. 

12  NAVSTA  -  FCRT  SERVICES 

13  NAS  MFAMAF  -  AIR  CPS 
1^  NAS  MFAfAF  -  AIMD 

15  fAS  KIPAMAP  -  SUPPLY 

16  NTC  -  FECTFACEN 

17  NTC  -  SEPSCCLCCM 

18  J-ANFCWEP  -  ACTIVE  CLTY 

19  MANPCVER  -  CEPS  CF  ACOL 

20  MANFC^EP  -  CIVILIAN 


Demand 

1057143,00    PECUISITIQNS 

79AU5.00    MAN-DAYS 

28877!5,00  MAN-hOURS 

6A788T8.00  MAN-HOURS 

194661.00  MAN-HOURS 

1074960,00  MAN-HOURS 

609420.00  GALLCNS 
3751626.00  MILLION  BTUS 
289142.00  ADMISSIONS 
498157.00  CUTPATIENTS 
45684.00  PERSCNNEL  ACTICNS 
9866.00  FCRT  SERVICE  HCLPS 
109424.00  AIR  OPERATIONS 
1499711.00  MAN-HOURS 
572336.00  REQUISITIONS 
5530.00  STUDENTS 
4716.00  STUDENTS 
759C4.00  PERSCNNEL 
17856.00  PERSCNNEL 
282C9.00  PERSCNNEL 
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COMPUTER  ALGORITHM 


PPCGRAM  TO  CQ  A  SENSITIVITY  ANALYSIS  OF  THE  NAVY  PERSONNEL  * 

RE5EAPCH  S  CEVELOPEMENT  CENTER'S  REGIONAL  INPUT-OUTPUT  * 

NCCEL  BV  MAKING  PARAMETRIC  CHANGES  TC  EACH  A(I,  J)  IN  THE  ♦ 

TECHNICAL  COEFFICIENT  MATRIX  AND  SCLVING   THE  RESULTING  * 

SYSTEM  CF  LINE/R  SIMULTANEOUS  EQUATIONS.   FOP  KNOWLEDGE  * 

GF  THE  MODEL,  SEE  "A  REGIONAL  INPUT-OUTPUT  MOCEL  FOR  FORE-  * 

CASTING  SHCRE-EASEC  NAVY  WORKLCAO,"  NPROC  TR  78-32,  AUGUST  ♦ 

* 

FPCGPAMMEC  BY  LCOR  LEE  T.  WOMACK,  lit  USN,  422-88-31A7/1 UC  * 

THESIS  PRCJECT  ♦ 

* 

ARRAYS:  * 

A       -  TECHNICAL  COEFFICIENT  MATRIX.  ♦ 

AINV    -  A-INVERSE.  ♦ 

AP      -  PARAMETRIZEC  TECHNICAL  COEFFICIENT  MATRIX.  ♦ 

^LSQ  USED  AS  WORKING  ARRAY  FOR  SCKE  * 

CCMPUTATIONS.  * 

BASE    -  MCDEL  SOLUTION  WITH  CELTA  FOR  P ARAMETRIZATI CN  ♦ 

ECUAL  TO  ZERO.  ♦ 

D       -  OFIGINAL  DEMAND  VECTOR.  ♦ 

DEL     -  VECTOR  OF  DELTA  VALUES  FOP  PARAMETRIZATION .  ♦ 

OUTPUT  -  MATRIX  CF  OLTPUTS  AT  EACH  ITERATICN  OF  * 

PARAMTERIZATION.  ♦ 

SECT    -  VECTOR  OF  ACTIVITY  NAMES.'  ♦ 

UNITS   -  VECTOR  OF  UMTS  FCR  EACH  SECTOR.  * 

WKAREA  -  WCRKAREA  FOR  LINV2F  TO  COMPUTE  "I-A  INVERSE."  ♦ 

XP4Y    -  VECTOR  OF  RCW  VALUES  CF  OUTPUT.   USED  TO  ♦ 

COMPUTE  CONFIDENCE  INTERVALS.  * 

SIERCLTINES:  ♦ 

CFIVER  -  RUNS  THE  ENTIRE  PROGRAM.  ♦ 

PARA    -  THE  HEART  OF  THE  ALGORITHM;  DOES  ROW  t    COLUMN  ♦ 

PARAMETRIZATIONS.  * 

LINV2F  -  CCMPUTES  "I-A  INVERSE."  ♦ 

XPRINT  -  KFITES  THE  A,  AINV,  AND  D  MATRICES.  * 

BMATRX  -  CCMPUTES  THE  INPUT  MATRIX,  B.  * 

UPRINT  -  VRITES  MATRICES  OF  ITERATIVE  OUTPLTS.  ♦ 

CONF    -  CCMPUTES  CONFIDENCE  INTERVALS  AND  WRITES  ♦ 

STJSTISTICS.  ♦ 

XCCNF   -  SETS  UP  THE  OUTPUT  MATRICES  TO  hAVE  CONFI-  ♦ 

CENCE  INTERVALS  COMPUTED.  ♦ 

COMPR   -  DCES  COMPARISONS  AND  WRITES  RESULTS  OF  PAPA-  ♦ 

METRIC  CHANGES.  ♦ 

MUL     -  DCES  MATRIX  MULTIPLICATION  TO  CCMFUTE  OUT-  ♦ 

PLT  (X).  ♦ 

TABLE   -  PRINTS  THE  (I-A)  AND  D  MATRICES.  * 

VARIABLES:  * 

IPRINT  -  PRINT  CCNTRCL  PARAMETER.  ♦ 

MSIZE   -  MATRIX  SIZE  (SQUARE  MATRIX).  ♦ 

NF      -  OSRN  FOR  RE/O.  * 
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NS      -  CSRN  FOR  WRITE.  * 

IWCRK   -  CIM6NSIGN  CF  WQRKABEA  FOR  LINV2F,  ♦ 

t*<lZE2    -  MSIZE  TIMES  TWO.  >» 

ZL      -  ZERO  TOLERANCE  LEVEL.  * 

lER     -  EPRCR  CCUNT  USED  IN  LINV2F.  ♦ 

LIL  -  NLMBER  OF  PARAMETRIC  ITERATIONS  CESIRED.  ♦ 
ICGLNT  -  CCL'NTS  ROW  ANC  CCLLMN  PARAMETRIC  ITERATIONS.     * 

ICEL7A  -  DELTA  VALUE  FOR  A  PARAMETRIC  ITERATION.  ♦ 

N       -  NUMBER  OF  OBSERVATIONS.  ♦ 

lOCF    -  CEGREES  OF  FREEDOM.  * 

TCTAL   -  SIM  OF  X«S.  ♦ 

XE4R   , -  MEAN  OF  X«S.  ♦ 

VARX   '-  VARIANCE  OF  X«S  (UNBIASED  ESTIMATE).  * 

SOEV  -  STANDARD  DEVIATION  OF  X«S.  * 
T       -  T-STATISTIC  FOR  CCNFIDENCE  INTERVAL  COMPUTATION.* 

XL      -  LCViER  CONFICENCE  LIMIT.  ♦ 

XL      -  UPPER  CONFICENCE  LIMIT.  ♦ 

CL      -  CONFIDENCE  INTERVAL  LENGTH.  ♦ 

NM     -  NLMBER  OF  CESERVATIONS  MINUS  I.  ♦ 

IPl     -  I  PLUS  1.  ♦ 

HOLD    -  CL'HHYf  USED  IN  SORT  ROUTINE.  ♦ 

ISTART  -  START  PARAMETER  FCR  CLTPUT.  ♦ 

ISTOP   -  STOP  PARAMETER  FCR  OLTPUT.  ♦ 

NN      -  N  DIVIDED  BV  2.  * 

KN      -  N  CIVIDED  BV  2»  TIMES  2.  * 

NNPl    -  N  CIVIDED  BY  2,  PLUS  1.  ♦ 

IN      -  N  DIVIDED  BV  4.  ♦ 

INN     -  N  CIVIDED  BV  4,  TIVES  2.  ♦ 

INPl    -  N  CIVIDED  BV  4»  PLUS  1.  * 

CI      -  LCWER  QUARTILE.  ♦ 

Q2      -  NEOIAN.  * 

Q3      -  UPPER  QUARTILE.  ♦ 

Kl      -  N  DIVIDED  BY  2f  PLUS  N  DIVIDED  BY  4,  PLUS  1.     * 

K2      -  N  CIVIDED  BV  2»  PLUS  N  DIVIDED  BV  ^,  PLUS  2.     ♦ 

RANGE   -  RANGE  OF  X'S.  * 

CRANGE  -  UTERCUARTILE  RANGE.  * 

CCOS    -  QLARTILE  COEFFICIENT  GF  SKEWNESS.  * 

KftNCl   -  INTEGER  ROUND  OF  Kl.   10  PERCENTILE  INDEX  OF     * 

X'S.  * 

KRND2   -  INTEGER  ROUND  OF  K2.   90  PERCENTILE  INDEX  OF     * 

X'S.  ♦ 

PCGK    -  PERCENTILE  COEFFICIENT  OF  KURTOSIS.  ♦ 

* 

CIMENSION  A(25t  25)t  AINV(25,  25),  AP(25,  25)t  BAS?(25),  C(25), 

1  CEL(25}»  0UTPUT(i5,  50),  SECT(25»  24),  UNITS(25,  20), 

2  *»KAR£A(75C),  XPA\(25) 

P?AC  THE  INITIALIZING  PARAMETERS. 
READ  (4,  101)  IPRINT,  MSIZE,  NF,  NS 

IVCRK  »  (MSIZE  ♦*  2)  +  5  *  >'SI  ZE 
^'SIZE2  »  2  ♦  MHZE 

CALL  CRIVER  (A,  AINV,  AP,  BASE,  D,  DEL,  OUTPUT,  SECT,  UNITS, 
1  WKAREA,  XPAV,  IPFINT,  IWORH,  MSIZE,  fSIZE2,  NF,  NS) 

hPITE  <NS,  102) 

STCF 

ICl  FCPNAT  (  715  ) 

1C2  FORMAT  (•  •,  ///,  •  ALL  CCMFUTATICNS  NOW  COMPLETE,  WITH  SUMMARY*, 
1  •  DATA  LISTED  ABCVE.   HAVE  A  NICE  DAY.'  ) 

E^C 

SIEBCLTINE  CRIVER  (A,  AINV,  AP,  BASE,  D,  DEL,  OUTPUT,  SECT,  UNITS, 
1  UKAPEA,  XRAV,  IPPINT,  IWORK,  MSIZE,  MSIZE2,  NF,  NS) 
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1  cA5  c(rSIZc)« 

2  SECTCVSIZEf    2 

ZL    a    l.OE-6 


JIZEj    MSIZE),    /iINV(MS!ZE,    MSIZE),    fiPiHSlZEi    HSMl 
C(MSIZEJf    OEL(MSIZE)f    Ol'TPUK  MSiZE ,    MSIZE2)i 
!^Jt    UNITS<MSIZE,    20)t    WKAREA( IWORK I ,    XRAY(MSI2E) 


fE>iC    IN    TJ-E    A    COEFFICIENT    MATRIX. 
CC    201    J    *    If     f'SIZE 
«      ^^    PE^O    (NFt    204)       (Ad,    J),    I   «    1,    MSIZE) 
^01    CONTINUE 

PEAC    IN    THE    OEKANC    MATRIX. 

READ    (NF,    204i       (C(I),    I    «lt    MSIZE) 

*  DI^'ENSICNALITY    OF    2A    COINCIDES    WITH   THE    DIMENSION         ♦ 

♦  GF    24    IN    SECT(MSIZE,    24).  * 

RE^C    EACH    SECTCP'S    IDENTIFICATION. 
CC    2C2    I    *    1,    fSIZE 

RE4C    (NF,    205)       (SECTd,    J),    J    «    1,    24) 
202    CCNTINUE 

^^tC    EACH    SECTOR' <    CNITS. 
CC    203    I    =    1,    NSIZE 

REieC    (NF,    206)         (UNITSd,    J),    J   «    1,    MSIZE) 

C^LL  BMATRX  (A,  AF,  0,  SECT,  ZL,  NSIZE,  NF,  NS) 

WFITE  THE  HE/iCING  AND  PRINT  THE  A  4ND  D  MATRICES. 

IF  (  IPRINT  .EC.  C  )    WRITE  (NS,  207) 

IF  (  IFPINT  .EC.  1  )    WRITE  (NS,  208) 

CALL  SPRINT  (A,  C,  SECT,  IPRINT,  NSIZE,  NF,  NS) 

C4LL  FARA  (A,  /IINV,  AP ,     EASE,  D,  CEL,  OUTPUT,  SECT,  UNITS, 
1  VKAREA,  XRAV,  ZL,  IPRINT,  IWORK,  MSIZE,  MSIZE2,  NF,  NS ) 

FETUFN 

204  FCPM^T  (7F1C.C1 

2C5  FCPN/T  (24A1I 

206  FCFf^T  (20A1) 

2C7  FCPNAT  ('ITHE  CRICIN^L  "A  AND  D"  MATRICES:*) 

208  FCRMAT  (UTHE  ORIGINAL  "A"  MATRIX:*) 

ENC 

SIEPCLTINE    PAR^     (/,    AINV,    AF,    BASE,    D,    DEL,    OUTPUT,    SECT,    UMTS, 
1    UK4REA,    XRAY,    ZL,    IPRINT,     IWORK,    MSIZE,    MSIZE2,    NF,    NS) 

CIMENSICN    A(MSIZEt    MSIZE),    4INV(MSIZE,    MSIZE),     iJP(MSIZE,    MSIZE), 

1  EASE(MSIZE),  C(^SIZE),  DEL(MSIZE2),  OUTPUT (MS IZE ,  MSIZE2), 

2  <ECT(NSIZ£,  24),  UNITS(MSIZE,  20),  UKAREAdWCRK )  ,  XRAY(MSI2E) 

CCMFLTE  THE  BASE  VECTOR  FOR  COMPARISON. 
DC  302  I  =  1,  VSIZE 

DO  301  J  *  1,  f-SIZE 
APd,  J)  =  4(1,  J) 

IF  (  I  .NE.  J  )    AP(  1,  J)  «  -APd,  J) 
IF  (  I  .EC.  J  )    APd,  J)  «  1.0  -  APd,  J) 
301     CCNTINLE 

2C2  CCNTINUE 

cm    dNV2F    ifiPi    ^SIZE,    MSIZE,    AINV,     l,    WKAREA,     IEP) 

CC    304    I    =    1,    HSIZE 
BASEd)    *    CO 

DC    203    J    =    1,    NSIZE 

IF    (    A8S(AINVd,    J))     .LE.    ZL    )         GO   TC    302 
BASEd)    a    BASE(I)     ♦   AINVd,    J)    *   0(J) 


98 


2C3     CCNTINUE 

3CA  CCNTIME 

PEAC  THE  NUMEEP  OF  CELTA  POINTS  ANC  THE  DELTA  VALUES. 

PEAC  (NF,  320)  LIL 

PEAC  <NF,  321)     (DEL(I),  I  «  1,  LIL) 

PRINT  TAELES  FCR  THE  THESIS  DRAFT. 

CALL  TABLE  <AF,  Ct  SECTt  UNITS,  MSIZE,  NF,  NS) 

CC  RCW  AND  COLUMN  PARAMETRICS  FOR  EACH  DELTA  VALUE. 
CC  31S  LEE  *  1,  LIL 

hRITE  (7,  222)  CEL(LEE) 

ICCLNT  =  0 

PARAMETRIC  ROW 

DC  210  I  a  1,  NSIZE 

CO  3C6  IL  «  1,  MSIZE 

CO  3C5  JL  «  1,  MSIZE 

AP(IL,  JL)  «  A<  IL,  JL) 

CONTINUE 

CONTINUE 

CG  207  J  :«  1,  MSIZE 

IF  (  ABS(AP(I,  J))  .LE.  ZL  )    GO  TO  307 
AP(I,  J)  «  AP(I,  J)  >  OEL(LEE)  *  AP(I,  Jl 

CONTINUE 

CO  2C9  11  «  1,  MSIZE 

CO  30e  JJ  «  1,  MSIZE 

IF  (II  .NE.  JJ  )    APdl,  JJ)  «  -APdl,  JJ) 

IF    (     II    .EQ.    JJ    )  APdl,    JJ)    »    l.C    -    APdl,    JJ) 

CCNTINUE 

CONTINUE 

CALL  LINV2F  (AP,  MSIZE,  MSIZE,  AINV,  1,  WKAPEA,  lER) 

CALL  fUL    (AINV,  0,  OUTPUT,  ZL,  ICOUNT,  MSIZE,  MSIZE2) 
CONTINUE 

FARAMTERIC  COLUMN 

CC  316  J  =  1,  KSIZE 

CO  212  IL  a  1,  MSIZE 

CC  311  JL  =  1,  MSIZE 

AF(IL,  JL)  =  A(  IL,  JL) 

CONTINUE 

CONTINUE 

CO  313  I  »  1,  MSIZE 

IF  (  AeS(AP(I,  J))  .LE.  ZL  )    GO  TO  312 
AP(I,  J)  3  AF(I,  J)  ♦  CEL(LEE)  *  AP(I,  J) 

CONTINUE 


99 


OC    315    II    »    1,    MSIZE 

00   314    JJ    «    1,    MSIZE  . 

IF    (     II    .NE.    JJ    )         AP(!I,    JJ)    a    -APdl,    JJ) 
^^    IF    (     II    ,£0.    JJ    )         /iPdl,    JJl    »    l.G    •    APdl,    JJ) 
314  CCNTINUE 

215  CONTINUE 

CALL    LUV2F    <AP,    MSIZEt    MSIZE,    AINV,    1,    WKAF6A,    lER) 

CALL    KUL    (AINV,    0,    OUTPUT,     ZL,     ICOUNT,    MSIZE,    MSIZE2) 
3U  CCMINLE 

IF    (    DEL(LEE)     .LT,    ZL    )         IDELTA   -    IFI X( (DEL (LEE)    *    lOC.C) 

IF  I  CEL(LEE)  •GT,  ZL  )    IDELTA  «  I FIX( (DEL ( LE E )  *  100. C) 
1    ♦  C.5) 

CALL  CCHPR  (iJP,  EASE,  OUTPUT,  SECT,  IDELTA,  MSIZE,  MSIZE2,  NF, 
1    NS) 

IF  i  IPRINT  .NE,  0  .AND.  IPRINT  .NE.  I  )    GO  TO  318 

WRITE  (NS,  223)  MSIZE,  MSIZE2,  IDELTA 

CALL  UPRINT  (OITPUT,  SECT,  ICOLNT,  MSIZE,  NF,  ^S,  1) 

317  ISTART  a  1 

IF  (  IPRINT  .EC.  2  )    GC  TO  31^ 

CALL  XCCNF  (QUIPLT,  SECT,  XRAY,  IDELTA,  ISTART,  MSIZE,  MSIZE, 
1    MSIZE2,  NF,  N<J 

ISTART  s  MSIZE  •»  1 

CALL  XCCNF  (OUTPUT,  SECT,  XRAY,  IDELTA,  ISTART,  MSIZE2,  f'SIZE, 
1    ^5I2E2,  NF,  NS) 

GO  TO  31«; 

318  IF  <  LEE  .EC.  LIL  J    GG  TO  317 
21S  CCNTINLE 

RETURN 

320  FCRMAT  (15) 

321  FCFM^T  (7F10.0) 

222  FCRMAT  (•  PARAPETPICS  WITH  A  DELTA  VALUE  OF',  F6.2) 
323  FCPMAl  CITHE  •,  12,  •  BY  «,  12,  •  PARAMETRIZED  RCVt«, 

1  •  ANC  COLUMN  fATPIX  V.ITH  DELTA  «•,  14,  •  PERCENT.',  ///) 

ENC 
UBPCLTINE  XPRINT  (A,  D,  SECT,  IPPINT,  MSIZE,  NF,  NS) 

CIMENSION  A(MSIZE,  MSIZE),  C(MSIZE),  SECT(MSIZE,  24) 

SLFPCLTINE  TO  WRITE  CUT  INTERMEDIATE  RESULTS  CF  "THESIS," 
WHICH  IS  A    PROGRAM  TO  CHECK  SENSITIVITY  IN  THE  I/C  MODEL. 

♦  IF  IPRINT  =  0,  PRINT  BOTH  THE  A  AND  0  MATRICES.      * 

♦  IF  IPRINT  a  1,  PRINT  CNLY  THE  A  MATRIX.  ♦ 

♦  IF  IPRINT  »  2,  PRINT  CNLY  THE  D  MATRIX.  * 

IF  (  IFPINT  .NE.  C  .AND.  IPRINT  .NE.  1  )    GO  TO  403 

SETLP  TO  PRINT  OLT  THE  MATRIX  PASSED  BY  COLUMNS  CF  10  ANC 
PCWS  CF  MSIZE,  UNTIL  THE  MATRIX  ll   COMPLETELY  PRINTED  WITH 
CCLUMN  HEADERS. 
I5TAP1  -  1 
ISTCP  *  IC 
4C1  IF  (  ISTOP  .GT.  MSIZE  )   ISTOP  »  MSIZE 


VFITE  (NS,  4061   (I,  I  «  ISTART,  ISTOP) 

CC  402  I 

WRITE 

4C2  CCNTINLE 


CC  402  I  »  1,  MSI2E 

WRITE  (NS,  4C7)    (Ad,  J),  J  «  ISTART,  ISTCF) 


/ 
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IF  (  ISTCF  .EC,  M3IZE)   GC  TO  403 
ISTAPT  =  ISTAPl  ■«.  10 
IMCF  »  ISTQP  4  IC 
WFITE  (NSf  4061 
GC  TC  401 

4C3  If  (  IPRINT  .NE.  C  .AND.  IFFINT  .NE.  2  )     RETURN 

ViFITE  CCT  THE  CEMAND  f'ATRIX  (MSIZE  X  1). 
4C4  UFITE  <NS,  40S) 

CC  4C5  I  =:  1,  NSI2E 

k>PITE  <NS,  410)   (SECTd,  J),  J  *  1,  24),  0(1) 

405    CCNTINUE 

FETliPN 

4C6  FCPMAT  (•  •,  ////,  8X,  12,  9U3) 

4C7  FCRMAT  (  'O*  ,  1CF12.5) 

4Cfi  FCPf-Al  (•!•) 

4C9  FCRMAT  <UTHE  CRIGINAL  DEMAND  VECTOR:',  //) 

41C  FCPr'AT  <«G«,24A1,  F16.2) 

END 

SIERCITINE  BMATRX  (A,  B,  D,  SECT,  ZL,  MSIZE,  NF,  NS) 

CIMcNSIGN  A(MSIZE,  MSIZE),  E(MSIZE,  MSIZ6),  0(fSI2E),  SECT(NSIZE, 
1  24) 

CC  5C2  I  a  1,  HSI26 

DC  501  J  =  1,  VSIZE 

e(  I,  J)  »  CC 

IF  I    AeS(A(I.  J))  .LT.  ZL  )    GO  TO  501 
8(1,  J)  '    A({,  J)  *  D(J) 
5C1     CCNTIME 

5C2  CCNTIME 

VFITE  (NS,  503) 

CALL  LPPINT  (6,  SECT,  MSIZE,  MSIZE,  NF,  NS,  3) 

PETLPN 

5C3    FCPMAT    CITHE    INFLT    MATRIX,    E:«,    //) 

ENC 

SIBRCLTINE  UFFINT  (ARRAY,  SECT,  NCCLS,  MROWS,  NF,  NS,  ITYPE) 

CIMENHCN  ARRAY(MPGWS,  NCOLS),  SECT(MROWS,  24) 

IF  (  ITYFE  .NE.  1  )    GC  TO  2 
ISTAPT  =  1 
ISTCF  =»  5 
ICCIN  =  5 
GC  TC  eci 

2  ISTJiPT  =  1 
ISTCP  =  10 
ICCLN  =  10 

6C1  IF  (  ISTCF  .GT.  NCCLS  )    ISTOP  «  NCCLS 

IF  (  ITYPE  .EC.  1  )  WRITE  (NS,  604)  (I,  I  =  ISTART,  ISTCF) 
IF  (  ITYPE  .NE.  1  )  WRITE  ( NS ,  610)  (I,  I  =  ISTART,  ISTQP) 
V»PITE  (NS,  605) 

CC  6C;  I  =  1,  NRCVS 

IF  (  ITYPE  -EQ.  1  )    WRITE  (NS,  606)   (SECTd,  J),  J  =  1,  24), 

1    (AFRAYd,  J),  J  =  ISTAPT,  ISTGP) 

IF  (  ITYPE  .EG.  2  )    WRITE  (NS,  608)   (SECTd,  J),  J  =  1,  24), 
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1    (/JPRAYd,  Jl,  J  s  ISTAP7,  ISTQP) 

IF  (  I7YPE  •EQ,  3  )    WRITE  (NSt  609)   (SECTd,  J),  J  =  1,  24) 

1  (APR^YCI,  J),  J  »  ISTART,  ISTOP) 
eC2  CCMIME 

IF  (  ISTOP  .EC.  NCCLS  )    GC  TO  6C3 
ISTART  s  ISTART  4  ICCUN 
ISTCF  =  ISTQP  4  ICQUN 
^PITE  (NSf  607) 
GC  TC  601 

6C2  FETIPN 

604  FCPMAT  (•  ••  4CX»  I2f  4120) 

6C5  FCP^AT  (•0«) 

606  FORMAT  CC,  24A1,  5F20.2) 

6C7  FCPNM  CIM 

608  FCRNAT  (•0»,  24A1,  10F10.2) 

6CS  FCPVAT  CON  24A1,  lOFlO.O) 

61C  FCPMAT  (•  •,  BOX,  12,  9110) 

E^C 

UePGCTINE    CONF    Ot    N,    NF ,    NS) 

CI^ENSICN    X(N) 

SIBPCCTINE  TC  CO^FllTE  A  951  CONFICENCE  INTERVALt  GIVING  THE 
LCV^ER  AND  UPPER  LIMITS,  FOR  THE  MEAN  CF  A  STRING  CF  NUMBERS. 

IF  (  N  .LE.  1  )    RETURN 
ICCF  »  N  -  1 
TCTAL  =»  0.0 
ZL  «  l.OE-6 

CC  701  I  =  1,  N 

TCT^L  «  TOTAL  ■*    X(I) 
701  CCNTINL'E 

XEAR  a  TCTAL  /  FLCAT(NI 
VAF>  »  0.0 

CC  702  I  «  1,  N 

VAPX  a  VAR>  ♦  <  X(I)  -  XEAR  )  ♦*  2 
7C2  CCNTIMJE 

VAR>  «  VARX  /  FLCndCOF) 
SCEV  »  VARX  ♦*  0.5 

IF  (  ICQF  .GT.   0  .AND.  ICOF  .LE.  3  )  T  «  3.H 

IF  <  ICOF  .GT.   3  .ANC.  lOOF  .LE.  7  )  T  »  2.36 

IF  (  ICOF  .GT.   7  .AND.  lOOF  .LE.  15  )  T  »  2.13 

IF  (  lOGF  .GT.  15  .ANC.  lOOF  .LE.  20  )  T  »  2.09 

IF  (  ICOF  .GT.  20  .AND.  IDOF  .LE.  25  T  »  2. 06 

IF  (  ICOF  .GT.  25  .ANC.  lOOF  .LE.  30  )  T  «  2.C< 

IF  <  lOCf  .GT.  30  .AND.  IDOF  .LE.  40  )  J    '  2. 02 

IF  (  ICCF  .GT.  40  )  T  «  2.CC 

XL  «  )iBAP  -  {  T  *  SDEV  /  (FLOAT(N)  ♦*  0.5  )) 
XL  a  >BAP  ♦  (  T  ♦  SDEV  /  (FLOAT(N)  ♦♦  0.5  )) 
CL  =»  >U  -  >L 

PCLTINE  TC  SORT  THE  ARRAY  AND  PRINT  OUT  THE  LIST  CF  NUMBERS 
IN  THE  SCRTED  ARRAY. 

NM  =  N  -  1 

CC  7C4  I  =  1,  NMl 
IPl  «  I  4  1 
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CC  7C3  J  «  IPli  N 

IF  (  X(I)  ,16.  X(J)  )   GO  TC  703 
l-CLC  »  XU) 
X(  1)  «  X<J) 
X(J)  «  HCLC 
7C3     CCMINUE 

7C4  CCNTIME 

IF  (  N  .GT.  ICC  I    GO  TO  7C6 

WFITE  (NS,  713) 

I57^R1  =  1 
IS7GP  3  IC 
7C5  IF  (  ISTCP  .GT,  N  )   ISTOP  «  N 

*iPI7E  (NS,  71^)    (  X(n,  I  =  IST/IRT,  ISTOP  ) 
IF  (  ISTCP  .EC.  N  )  GO  TO  7C6 
IST;iP7  a  ISTART  4  10 
ISTCP  *  ISTOP  4  IC 
6C  TC  705 

WRITE  OLT  THE  STATISTICS  CCNPUTED  4E0V6. 

7C6  WFITE  (NS,  715)  N,  IDOF,  TCTALf  XEAR,  VARX,  SOEV,  Tt  XL,  >U,  CL 

CCMFUTE  CUARTUES  AND  CCEFFICIENTS  OF  SKEWNESS  ANC  KURTOSIS. 

IS  N  EVEN  CR  CCC? 

NN  s  N  /  2 

KN  3  NN  !»  2 

NNPl  «  NN  >  1 

IF  N  CDC  ERANCh  TC  708. 

IF  i  KN  .LT,  N  )    GO  TO  706 

:^     t     *     *     *     * 

'•  N  EVEN.  ♦ 

:#  4  1^  :»  «  4 

IN  =  NN  /  2 

INN  a  IN  *  2 

TNPl  «  IN  +  1 

IS  N/2  EVEN  OP  GDC?   BRANCH  TO  707  IF  ODD. 

IF  (  INN  .LT.  NN  )    GO  TO  707 

N  EVEh  AND  N/2  EVEN.   COMPUTE  QUAPTILES. 

Kl  a  NN  4  IN 

Ki  a  Kl  +  1 

CI  «  (  X(IN)  4  X(INPl)  )   /   2.0 

C2  «  (  X(NN)  4  X(NNPl)  )   /   2.0 

C3  «  (  >(K1)  •»  X(K2)    )   /   2.0 

GC  TO  710 

N  EVEN  ^NC  N/2  CCC.   COMPUTE  QUARTILES. 
7C7  K2  s  NN  ■»  INPl 
CI  a  XINPl) 

02  «  (  X(NN)  4  X(NNPl)  )   /   2.0 
Q2    '    >{K2) 


GC  TC  71C 


i»  4  i#  «  :(i  « 

*      N  ODD.  ♦ 

i»  4  4>  «  «  ic 
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708 


C 

c 

c 
c 


716 


I^  »  hN  /  2 

I^^  «  IN  ♦  2 

INPl  »  IN  ♦  1 

IF  N/2  GDDt  BPflNCh 

If  <  INN  .LT.  ^^  ) 


TC  709, 
GO  TO 


709 


N  CCC  ^NC  N/2  EVEN,   COMPUTE  QUARTILES, 


Kl 

CI 
Q2 
C2 

GC 


TC 


NN  4  IN 
Kl  ♦  1 
(  X(IN) 
XNNPl) 
(  X(K1) 
710 


♦  1 

-♦  X(INPl)  )   /   2.0 

4  X(K2I  )   /   2.0 


7C9 


N  CCC  iSNO  N/2  ceo. 
K2  =  NN  ♦  INPl  ■»  1 
CI  s  XINPl) 
C2  a  XNNPl) 
C2  «  X(K2) 


CL'ARTILE  COMPCTATICNS 
CCEFFICIENTS. 


COMPUTE  QUARTTLES 


COMPLETE,  CCNPUTE  SKEWN5SS  £ND  KURTOSTS 


71C 


711 


712 


P/!hGE    s    X(N)    -    X(l) 
GRANGE    =    C2   -    CI 
IF    (    ieSiCP^NGE)     .LT.    ZL 
CCCS    «    CC3    ♦    <-2.C    *   02) 
IF    (     /ES(QRANGE)     .LT.    ZL 
KFNCl    a    (0.1    *    FLC^T(N)  ) 
KRNC2    '    (0.9    ♦    FLCAT(N)  ) 
IF    (    4eS(CP^NGE»     .LT.    ZL 
FCCK    3    (QR/iNGE    ♦    C.5)    /    ( 
IF    (    /8S(GRANGEI     .LT.    ZL 


)    GO  TC  711 

♦  Cll  /  GRANGE 
)    QCOS  3  0.0 

♦  C.5 

♦  C.5 

)         GO    TC    712 
X(KRND2I    -    X<KRND1)     ) 
)         PCOK    =    0.0 


VFITE  CLT  THE  CQMFLTED  STATISTICS  AND  STOP. 


UPITE 
1  FCCK 

PETLPN 


(NS,  716J   Xl),  X(M,  RANGE,  CI,  Q2,  Q3,  CRANGE,  QCCS, 


713 
71A 
715 


PCPM4 
FCPM^ 
FCFNA 
N 

CC 
TC 

VA 
SC 

•  T- 

•  LC 

•  LF 

•  CC 
FCP^4 

•  f*fi 

•  P4 

•  LC 

•  fE 

•  UP 
IN 


(•  THE  SCRTED  ARRAY  CF  NUNBERS:',  //  ) 

(•0«,  1CF13.2  I 

(•  S  ///,  •  STATISTICS:',  //, 

•-  17,  /, 


F  = 
TAL  =  • 
4N  -  «, 
FX  ==  •, 
EV  «  •- 


15,  /, 
,  F18.2, 
F19.2, 

F1S.2, 
F19.2 


/, 


STAT  =  •,  F17.2,  /, 

VER  CCNFlCENCE  LIMtT  = 


F20.2, 
/ 


PER  CCNFlCENCE  LIMIT  =  -, 

NFICENCE  INTERVAL  LENGTH  * 

T  (»OPININUM  s  • , 

XIMUf'  «  «.  F2C.2, 

^GE  «  •,  F22.2,  /, 

V»ER  CUARTILE  =  •, 

CIAN  '    • ,  F21.2,  / , 

FEP  CUARTILE  »  •,  F13.2,  /, 

lER-QUARTlLE  RANGE  »  S  F14.2, 

;PTILE  CCEFFICIENT 


F19.2,  /, 
F19.2,  //, 

/  , 


F13.2,  /, 


I 


OF  SKEWNESS 


a      • 


,  F9.4,  /, 


•  Cl\    .__  .-_  

6  •  PERCENTILE  CCEFFICIENT  CF  KURTGSIS  *  •,  F7.4,  ///) 

E^C 

SUBROUTINE    XCCNF    (ARRAY,    SECT,   XRAY,     TDELTA,    IST^PT,     ISTOP, 
1    KSIZE,    MSIZE2,    NF,    NS } 

C1^E^SIC^    APR4V(MSIZE,    MSIZE2),     SECT(ySIZE,    24),    XRAY(MSIZE) 

IF    (     ISTAPT    .LT.    ►SIZE    )         WRITE     (NS,    803)    lOELTA 
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IF    (     ISTART    .GT.    hSIZE    )         V»RITE    (NSt    804 1    IDELT/S 

CC    805    I    «    1,    NSIZH 
JJ    s    0 

CO    aCl    J    «    ISTART,    ISTCP 
JJ    a   JJ    4    1 

XR/V(JJ)    a    4RRAY( I,    J  J 
8C1  CCMINUE 

IF    (     I    .GT,    1    )         V^RITE     (NS»    805) 

fePITE    (NSf    806)         (SECTd,    JJ)»    JJ   =    1,    24) 

CALL    CCNF    (XPA\,    MSIZE,NF,    NS ) 

8C2  CONTINUE 

PETLPN 

aC3  FCRf<A7  CISTATISTICS  ON  PARAMETRIZED  ROWS  WITH  CELTA  »'t  14, 

1  •  P6FCENT.*,  ///) 
8C4  FCPMAT  ( • I  STAT  I  ST ICS  ON  PARAMETRIZED  COLUMNS  WITH  CELTA  »»,  14, 

1  •  PERCENT,',  ///) 
605  FCPNAT  CI') 
8C6  FCPMAT  («  *   *   ♦   DE>'AND  FCP •,  24A1 ,  •  ♦   *   *•,  ///) 

EhC 

SLERCLTINE  CCMFR  CAP,  BASE,  OUTPUT,  SECT,  IDELTA,  MSIZE, 
1  MSIZE2,  NF,  NS) 

CIMENSICN  AP(MSIZE,  NSIZE),  BASE(MSIZE),  QUTPUTCMSIZE,  MSIZE2), 
1  SECTiNSIZE,  24),  AVGROW(25),  AVGCCL{25),  PCTR(25),  PCTCC25) 

WRITE  (NS,  910)  ICELTA 

CC  9C3  I  s  1,  NSIZE 
TOTAL  «  CC 

CC  901  J  a  1,  MSIZE 

TCTAL  «  TOTAL  +  OUTPUTd,  J) 
9C1     CONTINUE 

A\GRQV.(I)  «  TCTAL  /  FLQAT(MSIZE) 
TCTAL  -  0.0 
ISTART  «  MSIZE  *  1 

CC  902  J  s  ISTART,  MSIZE2 

TCTAL  =  TOTAL  ♦  OUTPUTd,  JJ 
9C2     CONTINUE 

AVGCCLd)  '    TCTAL  /  FLOAT(MSIZE) 
9C3  CCNTINLE 

»»FITE    (NS,    911) 

CC    9C<    1=1,    MSIZE 

PCTR(I)    a    (     (AVGROWd)    /    BASE(I))    ♦    100.0    )    -    ICO.O 
fCTCd)    ■     (    (AVGCOL(I)    /    BASEd))    '*   100.0    )    -    100.0 
WRITE    (NS,    912)         (SECTd,    J),    J    «    1,    24),    BASE(I),    AVGRCWd), 
1         PCTP(I),    AVGCCKDt    PCTC(I) 
9C^    CCNTINLE 

WRITE  (NS,  913)  ICELTA 

CC  9C5  I  a  1,  NSIZe 
II  »  I  ■»•  MSIZE 

PCTRd)  «  (  (CLTPUTd,  I)  -  BASEd)  )  /  BASEd)  )  ♦  100. C 
PCTCCI)  »  (  (CUTPUTd,  II)  -  BASEd)  )  /  BASEd)  )  *  100.0 
WRITE  (NS,  914)   (SECTd,  J),  J  =  1,  24),  EASEd), 
1    CLTFUTd,  I),  FCTR(I),  OLTPUTd,  II),  PCTC(I) 
9C5  CCNTINLE 

CC  9C9  J  s  1,  MSIZE 
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JJ    a    J    ♦    MSIZE 

taPITE    (NS,    913  J    lOELTA 

WRITE    (NS,    S15J    Jf    Jf    J,    JJ 


CC    SC6    I    a    1,    NSIZE 

PCTRd)     »    (     (OUTPUTd,    J)    -    B/iSE(n     )     /    B/1SE(T)     )    *    ICO.O 
IF    (    lOELT/i    .EQ,    10    I  AF(I,    J)    =    PCTR(I) 

PCTCd)    «    (    (CUTPCTd,   JJ)    -   B^SEd)     I    /    E4SE(I)    )    *    3C0.0 
WRITE    (NS,    S14)  (SECTd,     K),    K    »    1,    24),    EASEd), 

1         ^^    GUTPUTd,    J),    PCTRd),    OUTPUTd,    JJ),    PCTCd) 
9C^  CC^TI^LE 

TCIP  »  0.0 
TCTC  «  C.C 


CC  907  K  =  1,  FSIZE 

TCTP  *  TCTP  4  PCTR(K) 
TOTC  «  TCTC  ♦  PCTC(K) 
FMEANR  «  TCTR  /  PLCAT(K) 
PMEANC  «  TCTC  /  FL04T(K) 
9C7     CChTINCE 

WFITE  (NS,  916)  PMEANR 

CC  9CE  K  «  1.  ►SIZE 

IF  <  ABS(PCTR(K))  .GT.  ABS(FMEANR)  )    WRITE  (NS,  917)  K, 
1       (SECT(K,  L),  L  »  1,  24),  PCTR(K),  PMEANR 
9CE     CChTIME 

909  CCKTINUE 

IF    (     ICELTA    .EC.    IC    »  WRITE    (NS,    918) 

IF    (     ICELTA    .EC.    10    )  CALL    UPRINT    (AP,    SECT,    MSIZE,    MSIZE, 

1    ^F,    hS,    2) 

FE^LP^ 

910  FCPM^T  dlPOW  4NC  COLUMN  MEANS  FOP  EACH  SECTOR  WITF  DELTA  «', 
1  14,  •  PERCENT.',  ///) 

911  FCRN/Jl  (•  •♦  'SECTOR*,  27X,  'BASE',  19X,  'ROW  MEANS  6X , 
1  •«  OIF  (RCWI«,  11>,  'col  MEAN',  6X,  '«  OIF  (COL)') 

912  FCFNM  dC,  2^A1,  5X,  F12.2,  12X,  F12.2,  7X,  Ff.2,  13X,  F12.2, 
1  7X,  F5.2) 

913  FCFHAT  CIEFFECT  CF  ROW  PARAMETRICS  ON  EACH  SECTCR  WITH', 

1  •  DELTA  «S  14,  '  PERCENT.',  ///,  «  SECTOR',  27X, 

2  'EASE',  lax,  'ROV.  CLTPLT',  5X ,  't  OIF  (SEC)S 

3  11>,  'CCL  OUTPUT',  5X,  '1  CIF  (CCL)') 

91^  FCPMAT  CO',  24A1,  5X,  F12.2,  12X,  F12.2,  7X,  F7.2,  IIX,  ^=12. 2, 
1  IX,  F7.2) 

915  FCPNAT  COROW',  13,  •  -  COLLMN',  13,  •  OF  OUTPUT', 

1  •  «  *  ♦  *  COL',  13,  •  -  CCLUMN',  13,  '  OF  OUTPUT.)') 

916  FCPMAT  CC,  //,  «  SUMMARY  FOR  PARAMETRICS  CN  PCWS:   MEAN  =  ', 
1  FE.4,  '  ?.') 

917  FCFM/T  CCPERCENTAGE  CHANGE  FCR  SECTOR  ',  12,  •,',  IX,  24A1,  /, 

1  •  CF  ',  F6.2,  '?', 

2  •  GREATER  IN  A6SCLUTE  MAGNITUDE  THAN  THE  MEAN  PERCENTAGE  CVANGE' , 
2  •  OF  ',  F6.2,  'I.') 

918  FCPMAl  CIMATPIX  CF  PERCENTAGE  CHANGES  WITH  DELTA  EQUAL  TC  ICIC, 
1  //i 

ENC 

SLBPCITINE    MUL    (MNV,    D,    OUTPUT,    2L,     ICOUNT,    MSIZE,    MSIZE2) 

CIMENHON    AINV(MS1ZE,    MSIZE),    D(MSIZ£),    OUTPUT  (MS  IZE ,    MSIZE2) 

ICCUNT    =    KGLNT    *    1 

CC    ICC2    I    =    1,    MSIZE 

aLTPUT(  I,     ICCLhT)    =    0.0 

CO    ICCl    J    »    1,    MSIZE 

IF    (    ABS(AINVd,    Jd     .LE.    ZL    )         GO    TO    1001 
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QUTPLTdt    ICCUNT)    *    OLTPUTd,    ICOUNT)     4    AINVd,    J)    *    C(J) 
1001  CCMINUE 

C 
10C2    CCMIME 


C 


FETLP^ 

E^o 

SIEPCLTINE    T/^ELE    (APt    D,    SECT,   UNITS,    MSIZE,    NF,    NS) 


CI^'£^SIC^  apcnsize,  msize),  o(msize),  sect(msize,  2^), 
1  lmts(?'Size,  20 

c 

C      k»FITE  TABLES  CF  "(I-A)"  ANC  "D«*  MATRICES. 

WPITE  <NS,  11C3) 
r 

CC  llCl  1  -  1,  MSIZE 

IF  (  I  .EQ.  11  .OR-  I  .EC.  21  )     WRITE  (NS,  1103) 
WRITE  (NS,  1104)  I,  (AP(I,  Jl,  J  »  1,  MSIZE) 

1101  CChTUUE 


I 


VRITE  (NS,  11C3) 


CC  11C2  I  =»  1,  MSIZE 
,  WRITE  <  NS  ,  11C5)  (SECTd,  J),  J  =  1,  24),  0(1), 

'      1    (UMTSd,  J),  J  a  1,  20) 

11C2  CCNTINCE 
C 

FETLP^ 
C 

1103  FCRMAT  ('1«,  /////) 

IICA  FCPMM  (•  'f  26X,  'ROWi*,  I2t  •)•♦  7(/,  4F14.5)) 

ilC5  FCPMAT  CCfAH,     24A1,  2X,  FIX. 2,  IX,  20A1) 
C 

E^c 
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